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ABSTRACT 


This  report  contains  computer  programs,  Instructions,  and  sample 
Input-output  data  for  antenna  pattern  synthesis  as  developed  In  the 
previous  report  "Computational  Methods  for  Antenna  Pattern  Synthesis." 
The  programs  are  valid  for  point  sources  arbitrarily  distributed  In  a 
plane,  and  for  pattern  synthesis  In  this  plane.  Included  are  programs 
for  synthesis  with  (1)  pattern  magnitude  and  phase  specified,  (2)  pat¬ 
tern  magnitude  only  specified,  (3)  these  two  cases  with  a  constraint 
on  the  source  norm,  and  (4)  the  first  two  cases  with  a  constraint  on  the 
source  quality  factor.  Also  Included  are  programs  to  compute  and  plot 
the  specified  and  synthesized  patterns. 
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1.  INTRODUCTION 


The  programs  used  to  compute  the  examples  of  Scientific  Report  No.  2, 
entitled  "Computational  Methods  for  Antenna  Pattern  Synthesis,"  are  described 
and  listed  In  this  report.  Each  program  is  accompanied  by  an  explanation  of 
the  Input  data,  a  verbal  flow  chart,  and  sample  Input-output  data.  Equations 
of  the  previous  report  (IJ  will  be  denoted  (2-  ),  and  equations  of  this  re¬ 

port  will  be  dei'Oted  by  a  single  number. 

The  programs  are  written  to  be  valid  for  N  point  sources  arbitrarily 
distributed  In  a  plane,  and  for  pattern  synthesis  In  this  plane.  However, 
for  the  examples  all  synthesized  patterns  are  radiated  by  10  point  sources 
equally  spaced  on  the  half  of  an  ellipse  on  one  side  of  the  minor  axis.  The 
major  axis  of  the  ellipse  is  twice  as  long  as  the  minor  axis.  The  first  and 
last  point  sources  are  at  the  ends  of  the  minor  axis.  Pattern  synthesis  is 
accomplished  hy  adjusting  the  complex  aii-'plltudes  of  the  10  point  sources  so 
that  the  pattern  radiated  bv  these  point  sources  is  as  close  as  possible  to 
a  specified  pattern.  'I'lie  amplitude  of  the  specified  pattern  Is  sketclic! 

below. 


Y 

1 

> 

4 

field  pattern 

:  X  V 

- S' 

; 

X  5 

r  X ' — sources 

Figure  1.  10  point  sources  and  specified  pattern. 


Tlj  J.  R.  Mauti  and  R.  F.  Harrington,  "Computational  Methods  for  Pattern 
Synthesis,"  Scientific  Report  No.  2  on  Contract  FI 9628-73-0-00*47 
between  Air  Force  Cambridge  Research  Laboratory  and  Syracuse  University. 
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II.  PATTERN  SYNTHESIS 


The  program  of  this  section  computes  the  vector  I  of  point  source 
strengths  given  by  (2-16) .  This  1  is  the  result  of  pattern  synthesis  with 
phase  specified.  The  program  prints  the  magnitudes  and  phases  of  the  ele¬ 
ments  of  ?.  The  program  also  punches  on  cards  the  real  and  imaginary  parts 
of  the  elements  of  ^  for  input  into  the  program  of  section  IX  of  this  report. 

Punched  card  data  Is  read  Into  the  program  according  to 

READ(1,30)  N,M,NK,NG 
30  F0RH\T(413) 

RE.VJ(1,33)(BK(I),  I-l,  NK) 

33  FORMAT (5E14. 7) 

READ  (1, 10)  (X  (I),  >1,  N) 

10  FORMAT(10F8.4) 

READ(1,10)(Y(I),  I-l,  N) 

READ(1,10)(PHI(I),  I-l,  M) 
no  35  JK-1,  NK 
DO  36  JG-1,  NG 
READ(1,10)(G(I),  I-l,  M) 

36  CONTINUE 
35  CONTINUE 

The  X  and  y  coordinates  of  N  point  sources  are  read  in  as  X  and  Y.  The 
specified  pittern  g^  of  (2-16)  Is  read  in  as  G  at  M  angles  PHI  (denoted 
by  ^  in  Figure  1)  In  degrees  In  the  plane  of  the  point  sources.  The 
index  JK  of  DO  loop  35  denotes  the  propagation  constant  BK(JK).  DO 

loop  36  makss  it  possible  to  deal  successively  with  NG  different  speci¬ 
fied  pattens  at  the  same  propagation  constant  BK(JK). 

Mlnlmjm  allocations  are  given  by 

COMPLEX  C(N*N) 
niMENSTON  LR(N) 

in  the  subroutine  LINEQ  and  by 
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COMPLEX  G(M),  T(N*M),  TC(N*M) 

TT(N*N),  TG(N),  F(N),  TF(M) 

DIMENSION  BK(NK),  X(N),  V(N),  PHI(M), 

TFM(M),  FM(N),  FP(N) 

In  the  main  program. 

Nested  DO  loops  14  and  15  store  the  matrix  element  given  by  (2-18) 
in  T((m-l)*N+n)  and  its  complex  conjugate  In  TC( (m-l)*N+n) .  Assuming 
that  [W]  Is  the  Identity  matrix,  nested  DO  loops  16  and  17  store  the 
matrix  (T*T]  of  (2-16)  by  columns  In  TT.  The  subroutine  LINEQ  inverts 
the  M  by  N  matrix  stored  In  TT.  Nested  DO  loops  19  and  20  store  lT*]gjj 
of  (2-16)  in  TG.  DO  loop  22  stores  f  of  (2-16)  in  F  and,  assuming  that 
[V]  appearing  In  (2-9)  Is  the  identity  ...at rlA,  ||f|P  In  FF.  Do  loop  22 
discovers  that  the  element  of  f  largest  in  magnitude  is  stored  In  F(J2). 

-> 

DO  loop  24  stores  the  complex  synthesized  pattern  [T]f  of  (2-17) 
in  TF  and  the  magnitude  of  the  synthesized  pattern  in  TFM.  l|[T]f|p  is 
accumulated  In  TFS,  the  pattern  synthesis  error  ||[T]f  -  In  E,  and 

llg^lp  In  CG.  Just  after  DO  loop  24,  O  of  (2-15)  is  stored  In  Q  and  E  is 
normalized  by  dividing  by  GG.  Statement  59  punches  F  on  cards  for  possible 
input  Into  the  program  of  section  IX.  DO  loop  41  normalizes  the 

element  of  f  largest  In  magnitude  to  unity  and  then  stores  the  magnitudes 
and  phases  in  degrees  of  f  in  FM  and  FP  respectively. 

The  propagation  constant  BK  of  the  input  data  Is  such  that  the  point 
sources  are  spaced  about  a  quarter  of  a  wavelength  apart  on  the  ellipse. 
Actually,  BK,  X,  and  Y  were  calculated  by  a  short  auxiliary  computer  program 
not  listed  in  this  report.  In  the  sample  input  data  the  magnitude  of  the 
specified  pattern  G  Is  given  In  Figure  1,  G  Is  real  and  G  Is  evaluated 
at  $  “  5°,  15°,  25°,... 355°  in  the  plane  of  the  ellipse.  The  numbers 
printed  out  under  the  headings  FM  and  FP  in  the  sample  output  were  used  to 
construct  the  first  two  columns  of  Table  1  of  Report  No.  2. 


9 


LIST  INN  l)H  H/MThKN  SYN1H(-SIS  PMIIRRAK 


//  1 0034, FI- ,  1  .41 ,  'MAllT?  t.lflF*  ,HEr,inN  =  ?00K 

//  FXFC  WATFIV 

//fUl.FTOZFOOl  no  SYS()in=W,nCK=<RECFM»F,RLKSI7F«Fn| 

//Gll.SYSIN  1)1)  » 

JijriH  MAUTZ  ,T  IMF  =  l  .FARFSeAO 

SUMRODTINF  L  I  NFO  (  LI.  .C  I 

COMPLEX  r,(  100) ,sti)R,sto,<;t,s 

DIMENSION  LRIbO) 
on  20  l=l,LL 
LR (  I  )  =  l 

20  CnNTINUb 
M1=0 

no  IH  M=1,LL 
K  =  M 

on  2  l=M,LL 
K1 eMl+I 
K2=M1+K 

IF ( ARS(REAL(C(Kl ) ) ) ♦ARS ( A IMARI C ( K 1 ) ) )-ARS( REAL  I C ( K?  )  I  )-ARS ( A IMAOIC 
1 (K2)  )  )  I  2,2,6 

6  K=I 

2  CONTINOF 
LS=LR<M) 

LR|M)=LR(K ) 

LK(K)=LS 
K2=M1+K 
STORsl ./C(K2) 

J1=0 

no  7  .1=1  ,LL 

Kl=Jl+K 

K2  =  .n+M 

STn=c (Ki) 

C(K1  )=r,(K2) 

C (K2  laSTO^STOR 
.11  =.I1+LL 

7  CONTINUE 
Kl =M1+M 

C  (M  '  =STOR 
no  11  I=1,LL 
IF( 1-M)  12,11,12 
12  K1=M1+I 
ST=C  IKI  ) 

C (Kl  laO. 

.11=0 

no  If  J=1,LL 
Kl=j; +1 
K2  =  .)l  +M 

C(K1  =C(Kl  )-C(K2)'«ST 
.11  =.ll  +LL 

10  CONT' NUF 

11  CONTINUE 
Ml =M- +LL 

18  CONTINUE 
Jl  =  0 

no  9  ,1=1  ,LL 
IE(J”LR(J))  14,8,14 
14  LR.I  =  l  R(.)) 

J2=(I.R.I-1  Ic-LL 

21  no  i:i  1  =  1, LL 
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KP=J2+I 

K1«J1+I 

S»C(K?) 

C(K2)=C<K1  ) 

C (Kl  )bS 
13  CONTINUE 

LR(J)=LR(LKJ) 

LR(LR.I)=LR.I 
!F(J-LR(J)I  lA.H.H 
fl  J1bJ1+LL 
N  CONTINUE 
RETURN 
END 

COMPLEX  U,U1 .CUNJO.GI 3ft  ) » T ( 360) , TC I  360 ) , TT 1 100 ) ,  TGI  10 ) , F ( 10 ) 
COMPLEX  TF(3ft) 

DIMENSION  RK( lOl.XIlOl.YI  10)  .PHI (36) .TFM(3ft),FM( 10) ,FP(10I 
U=<0..l.) 

P|B3.141b‘J3 

C=PI/lflO. 

RFAn(l,30)  N.M.NK.NG 
30  FORMAT (41  3) 

WHITE(3,3?)  N.M.NK.NG 

32  FORMAT! '0  N  M  NK  NG ' / 1 X  ,4  I  3  ) 

READ(l«33l(nK(I  ),|Bl,KiK) 

33  F0RMAT(5E14.7 ) 

WRITE(3.34) (BKII  ) ,  I  =1 ,NK) 

34  FCIRHATCOHK'/dX.BElft.?)) 

REAOd.lOl  (X(I  l,I=l.N) 

10  FOHMATdOF8.4) 

WR|TE(3,13)(Xd  ),I  =  1,N) 

13  FClRMAT('0X'/dX,10F«,4)) 

REAOd  ,101!  Yd  I ,  I  *1  .N) 

WRITE(3,37)  (Yd  l,I=l,N) 

37  FORMAT('0Y'/dX,10FH.4)  ) 

RFADd  dOKPHI  d  I  ,  I  =  1.M) 

WR ITE(3,11 ) (PHI d  I  .  I  =  1,M) 

11  FORMATCOPHI  '/(IX.lOFB.ft)) 

ZM  =  M 

DO  36  JK=1.NK 
HHbHK( JK ) 

J1=0 

DO  14  Krl ,M 
S1=PHI (K)wc 

r,s=cns(sn 

SN=SIN(S1 ) 

DO  16  1  =  1, N 
Ji=Jdl 

SI=RP#(X(I  )>1‘CS  +  Y(  I  l*SN) 

T(J1)=C0S(S1  )+U>i'SIN(Sl  ) 

TC(Jll=r,ONJG(T(Jl  )  ) 

16  CONTINUE 

14  CONTINUE 
.11=0 

DO  16  .1=1,  N 
DO  17  I=1,N 
J2  =  .) 

.13=1 
J1=J1+1 
TT( J1 )=n. 

DU  IH  K=1 ,M 
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TT(Jl)  =  TT(Jl)+TC(.t3)*T(J?) 

.I3  =  J3+N 
JZ=J2+N 

18  CONTINUE 
17  CONTINUE 
16  CONTINUE 

(ALL  LINEOIN.TTI 
Of)  36  jr,«l,Nf. 

LEAD! I  .10) IG( I ) . 1=1 
WK1TE(3,1?) (fid ),I=1,M) 

12  (OHMATI  '0fi'/(lX,10E8.4)) 
f'O  19  J  =  1,N 
.I2  =  J 

’fi( J)=n. 

DO  20  K=1 ,M 

‘  fi(  J)  =  Tfi(J)+TC(.(2)’»fi(K) 

.)2  =  .)2  +  N 
20  {.ONTINUE 

19  i;ONTINUE 
l=F=0. 

i2  =  0. 

DO  22  J=1.N 
d  J)=(). 

)1=J 

j(l  2  3  !  =1  ,M 

d  J)  =  F(.I)  +  TT(J1  )«‘TG(I  ) 

J1=J1+N 

23  ;nNT!NUE 
!;i=E(J)#f.ONJfi(E(J)  ) 

-E  =  FE  +  <51 

IF(S1.LT,S2)  fil)  TO  22 
52  =  S1 
J2  =  >l 

22  :ONTINUE 
.11=0 
TFS=n. 

F=0. 

fifi  =  0. 

DO  24  K=l ,M 
TF {K)=0. 

DO  25  1=1 .N 
.)1  =J1  +  1 

TF(K  )=TF(K  )  +  T(.ll  )*F(  I  ) 

25  CONTINUE 
S1=TF(Kl*C0N.lfi(TF(K)  ) 

TFM(K  )  =  .S(JKT  (SI  ) 

TFS=TFS+S1 

Ul=TF(K)-fi(K) 

E  =  E  +  U1*CUN.1G(U1  ) 

fifi  =  GG+G(K)*CUNJfi(G(K)  ) 

24  CONTINUE 
0=ZM*FF/TFS 
E=E/GG 

WRITE(3,26)  Gfi.FF.E.O 

26  FORMAT ( '0GG=' ,F14. 7, '  FF=',F14.7,'  E='.E14.7,'  0=*,E14.7) 
WRITE(3,?7) (F( I ),I=1,N) 

27  FORMAT! '0F'/( 1X,5F14.7)  ) 

59  WR ITE(2,3fl  )  (F( I  ),I  =  l,N) 

38  F0RMAT(5F14.7) 

01=1. /F(J2) 
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on  j  =  i  ,N 
F( 

FMK))=CABS(F(J)  ) 

S1=RFAL(F(J)) 

S2  =  AIMAri(  F(  .1  )  » 

FPt J)*ATAN?(S?.S1  ) 

FP(J)»FP(J)/C 

41  CONTINUE 

WRITE(3,4?MFM(  Jl  ,J=l,N) 

42  FORMATI 'OFM'/(  lX,10F7.3n 
WRITe(3,43)  (FP(  J)  iJM.N) 

43  FORMAT! 'OFP'/dX.lOFT.l  )  ) 

WRITE(3,2H)(TFin,I  =  l.M) 

2fl  FORMAT! 'OTF '/! IX, lOFfl. 41) 

WRITE!3,2'»)  !TFM!  n  .1  =1  ,M) 

24  FORMAT! '0TFM'/!lX,10Ffl. 4) ) 

34  CONTINUE 

35  CONTINUE 
STOP 
END 

.  .  SDATA 

10  34  1  1 

0.24nfl«82E+00 

0.0000  5.3408  10.4048  15.5027  14.2448  14.2448  15.5027  10.4048  5.3408  0.0000 

-10.0000  -4.4320  -8.4704  -6.2844  -2.5554  2.5559  4.2844  8.4704  9.4320  10.0000 

5.0000  15.0000  25.0000  35.0000  45.0000  55.0000  45.0000  75.0000  85.0000  45.0000 

105.0000115.0000125.0000135.0000145.0000155.0000145.0000175.0000185.0000145.0000 
205.0000215.0000225.0000235.0000245.0000255.0000245.0000275.0000285.0000245.0000 
305.0000315.0000325.0000335.0000345.0000355.0000 


5.0141 

0.0000 

3.8437 

0.0000 

2.3442 

0.0000 

1.7434 

0.0000 

1.4142 

0.0000 

1 .2208 

0.0000 

1.1034 

0.0000 

1 .0353 

0.0000 

1.0038 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

SSTOP 

/* 

// 


PRINTED  OUTPUT 

N  M  NK  NO 
10  34  1  1 

HK 

0.2408HR2F+00 


0.0000  5.3408  10.4048  15.5027  14.2448  14.2448  15.5027  10.4048  5.3408  0.0000 


-10.0000  -4.4320  -8.4704  -4.2844  -2.5559  2.5559  4.2894  8.4704  4.4320  10.0000 
PHI 

5.0000  15.0000  25.0000  35.0000  45.0000  55.0000  45.0000  75.0000  85.0000  95.0000 
105.0000115.0000125.0000135.0000145.0000155.0000145.0000175.0000185.0000195.0000 
205.0000215.0000225.0000235.0000245.0000255.0000245.0000275.0000285.0000295.0000 
305.0000315.0000325.0000335.0000345.0000355.0000 
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5.0191 

0.0000 

3.8437 

0.0000 

2.3442 

0.0000 

1.7434 

0.0000 

1.4142 

0.0000 

1 .220R 

0.0000 

1 .1034 

0.0000 

1 .0353 

0.0000 

1.003H 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000  0.0000 


f.G=  O.S^b^SO'Jt  +  O:'  FF=  0.  n3(S^?'iF+0?  E=  0.3121504F+00  0=  0.  1 25<J270E +02 
F 

-0.5421209E  00  0 . b224HV0F-01  0.11H7H10F  01  0.59ftqH31E  00-0. 14ft70<.flE  01 
-0.944R900E  00  0.6552457E  00  0.1432544F  01-0 . 34^3 1 7RF  00-0. 7954 1 R IF  00 
-0.9229400E-02  0.1055H90F  01 -0. 3044 194F  00-0. R 304R94F  00  0. RR3731RF-01 
0.1247218F  01-0.73H1020E  00-0.8539429F  00  0.5231934F  00  0 . 2333544E -0 1 

FM 


•  0.312 

0.742 

l.OOO  0.903  0. 

502  0.405  0.507 

0.717 

0.447 

0.300 

FP 

-38.3 

173.9 

0.0  -147.4  32.4  -122.3  37.1 

-124.8 

14.4 

149. R 

TF 

2.5445 

0.1147 

3.1540  -0.0913 

2.7784  -0.4051 

1.7283 

-0.4121 

0,7714 

0,0494 

0.479R 

0.5543 

0.449H  0.5719 

0.4R21  0.1754 

0.4518 

-0.1075 

0,3149 

-0,1359 

0,35R4 

-0.2250 

0.2700  -0.459R 

-0.0353  -0.5732 

-0.2900 

-0.4708 

-0.3143 

-0.3447 

-0.201R 

-0.4032 

-0.0808  -0.4H98 

-0.0044  -0.4700 

-0.0040 

-0,3379 

-0.1089 

-0.2402 

-0.2452 

-0.3117 

-0.2431  -0.5031 

-0.120R  -0.4144 

0.047R 

-0.5477 

0. 1494 

-0.4074 

0.2410 

-0.2431 

0.2951  -0.0445 

0.1449  0.1741 

-0.1479 

0.1827 

-0.4018 

-0.0090 

-0.5147 

-0.1244 

-0.4892  -0.0737 

-O.RR27  -0.0299 

-0.72R5 

-0.0479 

0.0543 

-0.0554 

1 .3298 

0.0454 

TFM 

2.5691 

3.1553 

2.8080  1.7747 

0.7732  0.7344 

0.R454 

0. 7044 

0.4444 

0.3449 

0.4233 

0.5333 

0.5743  0.5530 

0.4fl2R  0.4509 

0.4944 

0.4700 

0.33R0 

0.2437 

0.3944 

0.5477 

0.4262  0.549R 

0.4339  0.354R 

0.2987 

0.2294 

0.24fll 

0.4019 

0.5295 

0.4931 

0.8832  0.7317 

0.0792  1.3314 
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III.  MACNtTUDl-  PATTKRN  SYNTHKSIS 

The  pvonram  of  this  section  Implements  the  f  leld  miijinltude  pattern 
synthesle  described  In  section  IV  of  Report  No.  2, 

Punched  card  data  le  read  In  according  to 

RF.AD(1,30)  fI.M.NK.NG,N9 
30  FORMAT(5I3) 

RKAD(1,33)(BK([),  I-l,  NK) 

33  FORMATISFU.?) 

RFAD(1,10)(X(I),  1-1,  N) 

10  FORMAT (10F8. A) 

RFADd.lO)  (Y(l),  1  =  1,  N) 

RFAD(1,10)(PHI(I),  I-l,  M) 

DO  35  .IK-1  ,  NK 
DO  36  .10-1,  NG 
RKAD(1,10)(G(I),  1-1,  M) 

36  CONTINUF. 

35  CO.NTINUE 

The  X  and  y  coordinates  of  N  point  sources  are  read  In  .as  X  and  Y.  The 

JSm 

starting  value  of  h  e  Is  read  In  as  0  at  M  angles  PHI  In  degrees 

m 

In  the  plane  of  the  point  sources.  The  magnitude  of  G(m)  Is  the  magnitude 

th  J®m 

of  the  specified  pattern  at  the  m  angle.  The  angular  portion  e  of  G(m) 

gives  the  starting  value  for  B  .  The  Index  JK  of  DO  loop  35  denotes  the 

m 

propagation  constant  BK(JK)  DO  loop  36  makes  It  possible  to  deal  successively 
with  NG  different  specified  patterns  at  the  same  propagation  constant  RK(.1K). 
Step  2  of  the  Iteration  procedure  described  In  section  IV,  Report  No.  2,  Is 
repeatedly  executed  until  the  pattern  synthesis  erroi  t  of  (2-19)  ceases  to 
decrease,  but  not  more  than  a  maxlmimi  o'"  N9  times. 

Minimum  allocations  are  given  by 
CnMPI.F.X  (:(N*N) 

DIMENSION  l.R(N) 


( 
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In  the  liubroutlne  LINKQ  iitnl  by 


COHIM.I-X  C(M),  l'(N*M),  TC(N*M),  TT(N*N), 

K(N),  rF(M) 

DTMKN'SION  RK(MK),  X(N),  Y(N),  HHI(M),  TFM(M), 

11(M),  K(N9+l).  Q(N9+l),  FF(N9+1), 

FM(N),  Fr(N) 

In  tlie  naln  program. 

Nfsted  DO  loops  14  and  1 -i  store  the  matrix  element  given  by  (13-18) 

In  T((m-l)*N+n)  and  its  complex  eimiiie.ale  In  Tt;(iii- l)*’.'+n) .  As.sumlne  that 
[W]  Is  I  he  Identity  matrix,  nested  1)0  loops  16  and  17  store  the  matrix 
(f*T)  ol  (2-16)  hy  columns  In  TT.  Ihe  subroutine  l.INEQ  inverts  the  N  by  N 
matrix  stored  In  TT.  DO  loop  19  stores  (f*!)  ^t*  of  (2-16)  by  columns  in  TC. 

1)0  loop  40  .stores  h  of  (2-19)  In  H  and  accumulates  ihlF  in  GG.  Steps 
1,  2,  and  3  of  the  Iteriillon  procedure  described  In  section  IV,  Report  No.  2, 
are  exe<  uted  for  the  .19^*^  time  inside  DO  loop  41,  DO  loop  42  stores  i  of 
(2-16)  n  F  and  accumu late.s  il  f  In  FF(J2).  DO  loop  51  stores  [T)if  of  (2-19) 

In  TF,  (Tjfl  In  TFM,  !ilT]f|K  In  TFS,  and  the  pattern  synthesis  error  t  of 
(2-19)  n  SI.  DU  loop  51  also  changes  the  phase  ot  G  to  that  of  (T)f  .  Since 
J0m 

uppenrs  only  in  e  ,  the  Iteration  procedure  described  in  section  IV, 

Report  Do.  2, can  be  carried  out  without  calculating  explicitly.  Just  after 

DO  loop  51,  0  of  (2-15)  Is  stored  In  f)(J2)  and  the  relative  pattern  synthesis 

error  — ^  In  E(.I2).  Execution  remains  In  DO  loop  41  only  as  long  as  the 
h'  '• 

relative  error  coiulnue.s  to  dei nase. 

‘h.atpment  59  punches  F  on  cards  lor  possible  Input  into  the  program  of 
section  IX.  hi)  loop  55  finds  that  r(.12)  Is  largest  In  magnitude. 

DO  loop  56  normalizes  the  element  of  F  largest  In  magnitude  to  unity  and  then 
stores  he  magnitudes  and  phases  in  degrees  of  K  In  FM  and  FH  respectively. 

The  magnitude  of  the  lust  element  of  f  Is  0.235  in  the  sample  output  of 
the  i".ig!iltude  pattern  synthesis  program  as  opposed  to  0.234  in  the  table  of 
seel  Ion  IV,  Report  No.  2.  This  slight  difference  is  due  to  the  fact  that  the 
data  fo’  the  table  of  section  IV,  Report  No,  2,  was  computed  with  FDRTG  Instead 
of  UATF  V.  The  synthesized  patterns  were  also  slightly  different.  For  instance, 
the  ele'.eiuh  element  of  TFM  was  0.0811  with  fORKI  as  opposed  to  0.0816  with  W.\TFIV. 
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LIST  INC.  OH  MAGNITUnt:  HATTFHN  SYNTHfSIS  PROGRAM 


//  (0034 ,FF ,?OS. 1 «A I , 'MA0T7 ♦ JOF' ,RFGinN=?OOK 

//  FXeC  WATFIV 

//GO.FTOZFOni  ni)  SYSnoT  =  H,r)r,H=(HFr,FM=F,HLKSI7F  =  RO) 

//GO.SYSKY  Of)  # 

iJOH  MAIJTZ  .TIMF  =  1  .PAGES=*0 

SOBROUI  INF  LlNFO(LL.r,) 

COMPLEX  C( 1001  .STOR .STO.ST.S 
OIMFNSION  LR (SO) 

00  20  I=ltLL 
LW( I )=I 

20  CONTINUE 
M1=0 

00  IH  M=1,LL 
K=M 

on  2  I=M,LL 

KlrMl+I 
K2=M1 +K 

lF(AMS(«FAL(C(Knn  +  AHS(AlMAG(C(Kl)))-ARS(RFAL(C(K?n  )-AHS(  AIMAG  (  C 
1 (K2) )  )  1  ?,?,h 
A  K  =  I 

2  CONTINUE 
LS»LH (M) 

LH(M)=LR(K) 

LW(K)=LS 
K2=M1 *K 
S  rnR  =  l ./C (K2  ) 

.11=0 

no  7  J=1,LL 
K1=J1 ♦K 
K2  =  .ll+M 
ST0=C{K1) 

C(K1)=C (K21 
C(K2)=ST0*STUR 
.I1=.I1+LL 
7  CONTINUE 
K  1=Ml+M 
C(K1 )=STaH 
''0  11  I=liLL 
IF(I-M)  12.11.1? 

12  K1=M1*I 
ST=C(K1 ) 

C(K1 )=0. 

.11=0 

00  10  .1  =  1  .LL 

K1=.IH-I 

K2=J1+M 

C(K1)=C(K1 I-C{K2)*ST 
J1  =  .I1+LL 

10  CONTINUE 

11  CONTINUE 
M1=M1+Ll 

IR  CONTINUE 
J1=0 

no  q  .l  =  l,LL 
IF(J-LR(JM  lA.H.lA 
14  LR.I  =  LR  (.1) 

J2=(LR.I-nt‘LL 

21  no  13  1=1. LL 
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K2=J?+I 
Kl=.ll  +  I 
S=C (K?) 

c  (K2  )=r.(Kn 
C(Kl  i=s 
13  CriNT  INUI- 

LM ( JI=LRtLHJ) 

I  H  (LRJ  )=LR.I 
1F(J-LR(J1)  lA.H.H 
H  >i1=J1+LL 
4  ClINTINUh 
^  FTURN 
[  ND 

(.tlMPLEX  (J,U1  ,C0N.I(;,f.(3(S)  tT(  3ftO)  .TC  13^0)  ,TT(100I.F(10>.TF(3ft) 
ni  MEMSKIN  RK(10l,X(in),V(10(,PHl(3h).TFM(3(S),H(36)tF((Sl),0(An 
DIMENSION  FF (61)  .FMIIO) ,FP( 10) 
ll=(0..1.) 

I'l  =  3. 141643 
t.  =  PI/lRO. 

|!EA0(1,30)  N.M.NK  .NG.N*) 

30  f:flRMAT(6I3) 

NM1TE(3,3?)  N.M.NK  ,Nr.,N4 

32  ((IRMAn'O  N  M  NK  MG  N^'/lx.^n) 

IIFADI  1  ,33)  (HK(  I  )  ,  I  =1  ,NK) 

33  FORMA! (6E1 4.7  ) 

!RITE(3,34) (HK(  I  )  ,I=l,NK) 

34  l-ORMAK '0HK'/(1X,6E14.7)) 

)FAD( 1 ,10) (X ( I) . I =1 .M) 

10  OIRMATl 10FH.4) 

WRI TE(3,13)  (X(  I  )  ,  I  =1 ,N) 

13  •ORMAT('0X'/(lX,10Ffl.4)) 

3FAD( 1 ,10) (Y ( I ) , I =l,N) 

VRITE(3,37)  (Yd  ),I*1,N) 

37  =fmMAT('OY'/(lX,10FP.4)) 

!FAD( 1 ,10) (PHI ( I  )  , I =1 .M) 

HRITEO.l  I  )  (PHI  (  I  I  ,1=1  ,M) 

11  :()RMAT(  'OPHI '/(  lX.IOFH.A)  ) 

7M=M 

30  36  JK=l,NK 
BH  =  (3K(,IK  ) 

,11  =0 

DO  14  K  =  I  ,  M 
SI =PHI (K l^C 

cs=cns(si ) 

SN  =  SIN(S1  ) 
no  16  I = 1  ,  N 
,11  =J1+1 

SI =RR*(X ( I  )*CS  +  Y (  I  )*SN) 

i(  ji  )=r,os(si  )+ii*siN(si ) 

TC(.)1  )=r.l1N,lf,  (T  (,ll  )  ) 

16  OONTINOE 

14  CONTINUE 
.11  =0 

DO  16  ,1  =  1  ,N 
on  17  I  =  1  ,  N 
.12  =  .) 

.)3=  I 
.I1=.)141 
TT  (  .))  )  =  0. 
on  1 H  K=1 ,M 
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Tr(ji  )  =  iT(.ii)+Tr.(j^)#T(j?) 

>n=j3+N 
>I?  =  >I2+N 
IH  criNTiNiie 
1  7  CnNTINIIF 
IT.  CIINTINOF 

CALL  LINFO(N,TT) 

Jl*0 

.U  =  0 

nn  i«j  .1  =  1  .M 
nn  3<7  i  =  i,N 
.I3  =  J3*-1 
FI  I  )  =  TCI J3) 

3‘»  CONTINDF 
on  20  I=1.N 
.I1=J1  +  1 
TC  1 .11  )  =  n. 

.12=1 

no  3H  K=1.N 

TCl.ll  )  =  Tr,l.ll  l-fTT  (.I?)*FIK) 

.I2  =  .)2+N 
3H  CONI INUE 
20  CONTINIIF 
m  CONTINDF 

no  3T.  .i(;=i,Nr. 

R FAD  I  1  .10)101  I  I . 1  =  1  .M) 

WRITt(3,12)(0(l ),I=1,MI 
12  FORMAT  I '00'/ ( IX, inFB.4)  I 
00  =  0. 

no  <t()  .1  =  1  ,M 
S1=0I,II*C0NJ0(0(.))  I 
H(.|)=SORT  ISl  1 

oo=r,o+';i 

40  CriNT  INOE 

WRITE  I  3,47)  00 

47  FORMAT  ( 'ONORM  SOllARFO  OF  SPFCIFIFI)  PATTERN*  '  ,  F  1  4 . 7  ) 
FFI 1 )=0. 

F(1 )  =  1  . 

01 1 )=0. 

DO  41  .I<»  =  1  .NR 

.I2  =  JR  +  1 

FF(J2 )=0. 

no  42  J=1  ,N 

.)1=J 

F(J)=0. 

no  43  I =1 .M 

F(.l)=F(.))+TC(Jl  )*0(  I  ) 

.I1=J1  +  N 
43  CONTINUE 

FFIJ2  )=FF  (.)2  )+F(  .l)*CnNJO(F  (J)  1 
42  CONTINUE 
.11=0 
TFS=0. 

S1=0. 

no  51  .1  =  1,  M 
TFIJ )=0. 

DO  45  I =1 ,N 
J1=J1+1 

TF(.))  =  TF(.))+TI.)1  )*F(  I  ) 

45  CONTINUE 
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S?  =  rt-  LI  )  <T.ilN.I(K  I  HLI  I  ) 

Tf-M(.|)=S(JRI  (S?) 

I  H  I  1  HM(  J  I  .(rO.O.  I  (ill  TO  •>? 

H-S=  TRS  +  S? 

ri  .)  1  =  tH(,))/T(-M(  ,1  I  (♦TK  J  ) 
hy  S?  =  TfM(  J  )-HLI  I 
S I =S1 ♦S^OS? 

■il  fONTIMIIh 

0(j2)=;m*fhli?)/tfs 

F  1.12  )^S1 /fir, 

IF  iF(  121  .r.F  .K.mn  r.n  m  44 

4  1  C.IINTINIIF 

44  WR  I  TF  (  •i.'iR)  FF  (  I  I  ,  I  =1  ,  J2) 

‘v.  F(IKM41  ( 'ONORM  SOllAKFO  MF  F  '  /  (  1  X  t SF 14 . 7  )  > 

WW  I  TF  (  2.4hl  :F(  I  )  ,  I  =  1,.I2I 
4/,  FtlHMAl  (  'ORF  .ATIVF  FRKOR  '  /  (  1  X  ,  5F  14 . 7  )  ) 

WK  I  TF(  2,*)  3  1  10(1  1,1  =  1,. 12) 

*>3  FORMA  r  1  'OQ  •  /  (  1  X  ,  SF  14 .7  )  ) 

WR  I  TF  I  3,4H 1  (F (  I  I  .  I  =  1 ,N) 

4R  format  I  'OF' /  (  1  X,  “iFK./ )  ) 

•  S4  WR  I  TF  I  2 , 33) (F ( I)  . I  =  1 ,N) 

S2  =  0  . 

IK)  S')  )=l,N 

r,  1  =F  i.i  I  »f.(iN,ir, I F 1 ,1 ) ) 

IF  (SI  .L  1  .S2  1  (ill  Til  SS 
S2  =  S  1 
.12  =  .) 

SS  rONTIMIF 
m  =  1  ,/F  I J2) 
nn  S4  .i  =  i,N 
F(  .)  )  =  III»F(.)  ) 

FM(J)=f,AHSIF(.ll) 

51  =RFAL  (Fl.l  I  ) 

52  =  AlMAr.(F(.l  )  I 
FP(JI=(ATAN2(S2,S1 ) I/C 

S4  CONTINUF 

WR  I  TE  (  3,57  I  (FM(J)  ,.1  =  1  ,N| 

57  FORMAT! 'OFM' /( 1X,10F7.3I ) 

WRITF(3,5fll (FP( J|,J=1,N) 

5«  FORMAT (  'OFF '/( 1 X , 10F7, 1 )  ) 

WR  I  TF(  3,49|  (TF(.II,.)  =  1,M) 

49  FORMA  I  1  'OTF ' / ( IX , 10FH.4 )  I 
WR  ITF  (  3,50)  (TFMIJ  )  ,  .1  =  1  ,M) 

50  FORMAT! 'OTFM'/nX.lOFfl. 4)  ) 

3ft  CnOTINUE 

35  CONTINUE 
S  I  OP 
FNO 

tllAT  A 

10  3ft  1  1  40 

0.290HRR2'-+00 

0.0000  ‘j.3ftOH  in.ftOAM  15.5027  lP.?ft4R  I9.2ft4R  15.5027  10.ft04R  5.3ftOR  0.0000 

-10. 0000  -9.ft320  -fl. 470ft  -ft.2fl‘74  -2.5559  2.5559  ft.2R94  R. 470ft  9.ft320  10.0000 

5.0000  15.0000  25.0000  35.0000  45.0000  55.0000  ft5.0000  75.0000  H5.0000  95.0000 

105.0000115.00001 25.00001  35 . 0000145 .0000155 .00001 ft5.  OOOOl 75. 00001 fl5 . 00001 95 .0000 
205.000021  s.  000022 5. 000023 5. 0000245. 0000255. 00002ft 5. 000027 5. 00002R5. 0000295. 0000 
30ft .0000315.0000- 25.0000335.0000345.0000355.0000 

5.0191  0.0000  3.8ft37  0.0000  2.3ftft2  0.0000  1.7434  0.0000  1.4142  0.0000 

1.220R  0.0000  1.1034  0.0000  1.0353  0.0000  1.003fl  0.0000  0.0000  0.0000 

n.nonn  o.onno  o.oooo  o.oooo  0.0000  o.oooo  o.oooo  o.oooo  o.oooo  o.oooo 
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0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000  0.0000 
ISTOP 
/* 

// 

PklNlEO  OUTPUT 
N  M  NK  NC.  Nk 
10  36  1  1  40 

BK 

0.2<»OflBH2E  +  00 
X 

0.0000  ‘).360fl  10.604H  IS.BOP?  l‘J.?64«  1<J.264S  16.3027  10.604R  6.360R  0.0000 

Y 

-10.0000  -‘*.6320  -8.4706  -6.2844  -2.6664  2.6664  6.2R44  R.4706  4.6320  10.0000 

PHI 

6.0000  15.0000  26.0000  36.0000  46.0000  66.0000  65.0000  76.0000  85.0000  45.0000 
106. 0000 116. 00001 26.0000136.000014  6.0000166.0000166.0000176.0000186. 0000146. 0000 
205 .0000216.0000226.0000236.0000246.0000266.0000266 .0000276.0000286.0000245,0000 
306.0000316.0000325.0000336.0000346.0000366.0000 


r, 


6,0141 

0.0000 

3.8637 

0.0000 

2.3662 

0.0000 

1.7434 

0.0000 

1.4142 

0.0000 

1 .2208 

0.0000 

1.1034 

0.0000 

1.0353 

0.0000 

1.0038 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.oooo 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

NORM  SQUARED  OF  SPECIFIED  PATTERN=  0.5654604E+02 


NURM  SQUARED  D 
O.OOOOOOOE+00 
0.266R625E+02 
0.2R60146E+02 
0.2760438E+02 

RELATIVE  ERRUR 
O.lOOOOOOE+01 
0.223R368E+00 
0.1727168E+00 
0.1720466E+00 


O.OOOOOOOE+00 

0.2237406E+02 

0.2242101E+02 

0.2152713E+02 

F 

-0.61  31  346E  +  00' 


F 

O.I 336623E+02 
0.3068034E+02 
0.2745184F+02 
0.2751226F+02 


0.2R481H4f +00 
0.2103622F+00 
0.1722322F+00 
0.1720466F+00 


0.1254267F+02 

0.2638416E+n2 

0.21R8480F+02 

0.2153314E+02 


0.2752376F+00 


0.1R14644F+02 

0.3176707E+02 

0.2765762E+02 

0.2762267F*02 


0.2668424F+00 
0. 14537R2E+00 
0.17212R7E+00 
0. 172044HF+00 


0.1604077F+02 

0.25H2610F+02 

0.2164H62F+02 

0.2154136F+02 


0. 1244744F+01 


0.2036316F+02 

0.311655HF+02 

0.27542RRE^02 

0.2753236E+02 


0.24471 lOF +00 
0. 181 7446E+00 
0.172in53E+00 
0.1720446F+00 


0.1768746F+02 

0.24B6603E+0P 

0.2166760F+02 

0.2154RR3F+02 


0.4462R00F+00 


0.2223022E+02 

0.2476671E+02 

0.2750768F+02 

0.2764047F+02 


0.2361352F+00 
0. 1 74R158F+00 
0.1720440F+00 
0. 1 720446F+00 


0. 1R46R51F+02 
0.2344670F+02 
0.2162474E+02 
0.215552RF+02 


0. 1632R41E+01 
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-0. 1 ‘)HHA70K  +  ni  O.-iYHT^lHfi+OO  0.  2502<S43P  +  01  0.  ?R  1 1  ?50e  +  00-0.  1  S'??  1  94E  +  0 1 
-0.3373H44E  +  00  0 . 1 3ft«»034F  +  01  0 .223#*3«  lF  +  00-0. 14355R2E  +  01  0,  lft07297F-01 
n.l  /imHnF  +  Ol-O.AlS'lYHOF+OO-0. 1279405F  +  01  0.4fll  2A00E+OO  0.35441 19E+00 

EM 


0.244 

0.424 

0.985  1.000  0.434  0.553  0.570 

0.444 

0.528  0.235 

EM 

125.0 

•r. 

• 

O 

1 

153.2  0.0  -159 

.2  24. 

7  -140.3 

10.3 

172.9  -42.8 

TE 

3.0114 

-0.0309 

3.4498  -0.4797 

2.5152 

-1.4827 

0.5754 

-1.4077  -1.1052 

-0  •  544f 

-1  .3524 

0.9  370 

-0.3444  1.5845 

0.5740 

1.0149 

0.4059 

0.1924  0.2243 

-0.0425 

0.077H 

0.0?47 

0.0547  -0.0043 

-0.1144 

-0.0908 

-0.2980 

-0.0544  -0.2747 

-0.022C 

-0.1144 

-0.1 273 

0.0100  -0.2474 

0.0312 

-0.1984 

-0.0444 

-0.0114  -0.1444 

0.0945 

-0  .2082 

-0.0335 

-0.0924  -0.1487 

0.0510 

-0.1782 

0.0181 

-0.0895  -0.0958 

-0.1372 

-0.0114 

-0.2322 

0.1820  -0.1449 

0.1350 

0.0754 

-0.0859 

0.0453  -0.1245 

-0.1401 

-0.0508 

-0.1402 

-0.2284  0.1054 

-0.5787 

0.2217 

-0.5413 

0.1203  0.2742 

0.050^ 

1.4839 

0.0144 

1EM 

3.0120 

3.5161 

2.9197  1.7074 

1.2420 

1.4453 

1.4220 

1.1487  0.4357 

0.2343 

0.0814 

0.0351 

0.1478  0.3029 

0.2775 

0.1725 

0.2474 

0.2011  0.0459 

0.  IPOf 

0.2082 

0.1  125 

0.1854  0.0913 

0.1474 

0.2525 

0.2339 

0.1547  0.0971 

0.202t- 

0.1491 

0.2317 

0.4197  0.5545 

0.2808 

1.4847 
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IV.  EIGENVALUKS  AND  F.TCKNVFICTOKS  OF  [T*T] 

The  program  of  tills  section  punches  on  cards  the  eigenvalues  X  and 
the  eigenvectors  of  (2-35)  when  W  and  V  are  Identity  matrices. 

Punched  card  data  is  read  In  according  to 

READ(1,30)  N,M,BK 
30  FORMAT (21 3,  E14.7) 

READ(1,10)  (X(I),  I-l,  N) 

10  FORMAT (10F8. 4) 

READ (1, 10)  (Y  (I),  1-1, N) 

RFJ1D(1,10)(PHI(I),  I-l,  M) 

The  X  and  y  coordinates  of  N  point  sources  are  read  In  as  X  and  Y.  The  M 

angles  of  (2-18)  are  read  In  as  PHI.  The  propagation  constant  Is  RK. 
n 

Minimum  allocations  are  given  by 

COMPLEX  A(N*(N+l)/2),  R(N*N) 

in  the  subroutine  CEIG  and  by 

COMPLEX  T(N*M),  TC(N*M),  A(N*(N+l)/2) , 

TT(N*N),  P(N*N) 

DIMENSION  X(N),  Y(N),  PHI(M),  AM(N) 

In  the  main  program. 

Nested  DO  loops  14  and  IS  store  the  matrix  element  given  by  (2-18) 

In  T((m-l)*N+n)  and  Its  complex  conjugate  In  TC((m-l)*N+n) .  Assuming  that 
[W]  is  the  Identity  matrix,  nested  DO  loops  16  and  17  store  the  upper  tri¬ 
angular  portion  of  [T*T]  of  (2-35)  in  A  and  the  whole  matrix  IT*T]  by 
columns  in  TT, 

The  subroutine  CEIG  called  by  statement  40  puts  the  eigenvalues  of 
[T*T]  In  decreasing  order  In  the  diagonal  positions  of  A  and  stores  the 
1^^  component  of  the  J^*’  eigenvector  of  [T*T]  In  P((j-1)*N+1) .  The  eigen¬ 
vectors  calculated  by  CEIG  satisfy  the  orthogonality  relationships  (2-36) 
and  (2-37).  If  the  fourth  argument  MV  of  CEIG  Is  not  1,  the  eigenvectors 
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are  stored  In  the  second  argument  R.  If  MV»1,  no  eigenvectors  are  stored. 
CEIG  Is  an  extension  to  Hermit  Ian  matrices  of  the  subroutine  EIGEN  In  the 
IBM  Scientific  Subroutine  Package,  f2],  {3],  A  detailed  description  of 
CEIG  Is  nut  Included  In  this  report. 

DO  ..oop  38  stores  the  eigenvalues  in  AM.  Statements  41  and  42 
punch  the  eigenvalues  and  eigenvectors  on  cards.  To  check  the  eigenvalues, 
DO  loop  2d  stores  In  AM(J)  the  quadratic  form  i|>j[T*T]^j  which  would  be  the 
Rayleigh  i|uotlent  for  1  j  if  ■  1. 

According  to  the  discussion  in  section  VII,  Report  No.  2,  the  eigen¬ 
vectors  111  the  sample  printed  output  are  supposed  to  be  purely  real.  How- 
ever,  becm&e  of  round  off  error,  the  computed  elements  of  [T*T]  and  thus 
the  eigenvectors  have  small  Imaginary  parts. 


[2]  IBM  5ystem/360  Scientific  Subroutine  Package  (360A-CM-03X)  Version  III, 
Programmer's  manual. 

[3]  C.  Frbberg,  "Introduction  to  Numerical  Analysis,  Addison-Wesley ,  1965, 
pp.  109-112. 
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LISTING  OF  PROGRAM  TO  PUNCH  RIGENVALUES  AND  EIGENVECTORS 

//  (0034,FFi I5S.lf42I t'MAUTTtJOE* .REGI0N-200K 

//  EXEC  WATFIV 

//Gn.FT02F001  DO  SYSOUT •RiDCR" I ReCFM»F tRLKS I / E»flO I 
//GO.SVSIN  DO  • 

$JOB  MAUTZ. TIME-1, PAGES-40 

SURROtJTINE  CEIG(A,R,N,MV) 

COMPLEX  A(55),R(100I 
COMPLEX  U,U1,U2,U3,C0NJG 
IJ=tO.,l  .) 

MVl-MV-1 

IF(MVl)  10,2S.10 
10  IJ«0 

DO  20  .1  =  1  ,N 
J1*IJ+J 
00  21  1=1, N 
I.I=IJ  +  1 
R  (  IJ  1=0. 

21  CONTINUE 
R(J1  »  =  1. 

20  CONTINUE 
25  AN0RM=0. 

IA  =  1 

DO  35  J*2,N 
Jl-J-1 

DO  36  1*1, J1 
1  A-IA  +  1 

ANORM-ANURM-fAl  I  A)«CONJG<  At  lAI  ) 

36  CONTINUE 
I A=1 A+1 
35  CONTINUE 
FN*N 

IF(ANORM)  165,165,40 
40  AN0RM*1.414*S0RT(AN0RM) 

ANRMX=AN0RM«l.E-6/FN 
IND=0 
THR*AN()RM 
45  THR=THR/FN 
127  LM-0 
51  IM0»0 

DO  55  M=2,N 

IMO-IMO+N 

J1*M-1 

LM-LM+1 

M0*LM 

MMsMO+M 

LL*0 

1L0=-N 

DO  50  L=1,J1 

ILO=ILO+N 

LM*LM+1 

LO=LL 

LL=LL+L 

RA«REAL I AtLM ) ) 

AA=AIMAG( AtLMI ) 

AM2  =  RA*RA'fAA*AA 
AM>S0RT( AM2) 

IF(AM-THR)  50,65,65 
65  IND*1 
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ALL. REAL! A(LL ) ) 

AMM.REAL(A(MM| ) 

X=.5*(ALL-AMM) 

SO=SORT( AM?+X*X ) 

SlNX=-SORT(AM?/(  SO+AHS<X ) ) ) ) 

IF(X(  70,7^1,75 
70  SINX.-SINX 
75  SINX2=SINXfSINX 
C0SX2  =  1  .-SINX2 
COSX=SOHT(COSX2I 
SINCS  =  SINX1<CC1SX 
SINTr*A/AM 
COST.RA/AM 

Ul=SINX*(r,OST+U*SINT> 

U2=C0NJG(U1 ) 

10  =  0 

on  12<1  1  =  1, N 

iF(i-L)  ao.n^.Hi 
HO  IF (  IMM)  ft2. 115,83 

82  IM=I+MO 
!L=!+LO 

U3=A(  IL  l»CriSX-A(  IM)*t)2 
A(  IM)  =  A(  IL  1*U1+A(  IM  KCnSX 
GO  TO  12A 

83  IM=M+IO 
IL=I +L0 

II3  =  A(  IL  )*(:OSX-cr)NjG(A(  IM)  )*II2 
A(  IM)=cnNJG(A(lL))*U2+A( IMI^COSX 
GO  TO  124 

Hi  IF(IMM)  84,115,85 

84  IMsI^MO 
IL.L+IO 

t)3  =  A(  IL  l*COSX-CONJG(A(  IMn*Ul 
A(  IM|=CONJGl  A(  IL  I  t*in+AI  1M»*C0SX 
GO  TO  124 

85  IM.M+IO 
IL=L+IO 

ll3=A(IL)*CnsX-A(IM)*Ul 
A(  IM)=AI IL)*U2  +  A( IM)*CnSX 

124  A( IL )=U3 

115  IF(MVl)  120,125,120 
120  ILR.ILO+I 
IMR=IMO+I 

U3=R ( ILR )*COSX-R  t IMR )*II2 

R  (  IMR  l=R  (  lUR  )*lll+R  ( IMR  )*CnSX 

R(ILR)=t)3 

125  I0=104l 
12ft  CONTINUE 

X=2.*AM*SINCS 

Y=ALL*C0SX2+AMM*SINX2-X 

X=ALL*SINX2+AMM*C0SX2+X 

A(LM)=A(LM)*C0SX2+II1'»(  (  ALL-AMM)*C0SX-C0NJG(A(LM)  )*Ml) 
AILL  )=Y 
A(MM)=X 
50  CONTINUE 
55  CONTINUE 

IFIIND-l)  160,155,160 
155  IND=0 

GO  TO  127 
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160  !F(THR-ANRMX)  165,165»45 
165  I0«-N 
LL-0 

on  lfl5  I»1.N 

JO-IO 

IO«IO>N 

MM*LL 

LL-LL+I 

DO  185  J-I,N 

JO"JO+N 

MM«MM+J 

ALL-REAL(A(LL) ) 

AMM>REAL(AIMH) ) 

IF(ALL-AMM)  170fl85.185 
170  A(LL)=AMM 
A( MMI=ALL 

IF(MVl)  175,185,175 
175  on  180  K«1,N 
ILR-IO+K 
IMRiJO+K 
U3cR(ILR I 
R(  ILR)>RI IMRI 
ISO  R ( IMR l«U3 
185  CONTINUE 
131  RETURN 
END 

COMPLEX  U,Ul,CONJ0,T(360»,TCI360)  ,A(55)  ,TT(100),P(100) 
DIMENSION  X(10),VI10),PHI  (36),AM(10> 

P I-3.1415R3 
C«PI/1R0. 

U=(0.,1.) 

REA0(1,30)  N,M,nK 
30  F0RMAT(2I3.E16.7| 

WRITEI3,32)  N,M,RK 

32  FORMATCO  N  M '  , 6X , 'BK • /IX  , 2 1  3 .F 14 .7 ) 

REA0I1,10)IX(I l,I°l,N| 

10  FORMAT! lOFR. 4) 

WRITE(3,13I(X(I ),I<1,N1 

13  FORMAT! ' OX • / ! 1 X , 10F8 .4 »» 

REA0!1,10I!Y!I ),I>l,N) 

WRITE!3,37MV!  I  |,I*l,Nl 

37  FORMAT! •0Y'/!lX,I0Ffl. 4)1 
REA0!1,I0)!PHI!I ),I>1,M) 

WRITE!3,11I!PHI !I I, 1*1, Ml 

11  FORMAT! 'OPHI •/! IX, 10F8. 41 1 
J1«0 

DO  14  K*1,M 
SI-PHI  IK 1*C 
CS-COS!  SI  I 
SN-SINISll 
DO  15  I-1,N 
Jl-Jl+1 

S1-BK*!X!I1»CS+Y!II*SN1 
T! Jll-COSISl I+U^SINISI 1 
TC! J1I=C0NJG!T!J1 ) ) 

15  CONTINUE 

14  CONTINUE 
Jl-0 
J5-0 

DO  16  J=1,N 


27 


.I4  =  J 

nn  17  i=i,j 

J2*J 

j^=i 

Jl=Jl>i 

A(  .11  1=0. 

on  ifl  K*i  ,M 

A(Jl)  =  A(,Jl)*Tr,(J3)-»T(J?) 

J3=J3+N 
J2=J2+N 
18  CONI  INtlE 
J8=JS+I 
TT(J6)=A( J1  ) 

TT(  J4)=CI)N.IG(A(  J1  1  ) 

J4=J4+N 
17  CONTINUE 
jS  =  Jb-fN 
liS  CONTINUE 

40  CALL  CEI  r.(  A.P.N.O) 

.11=0 

DO  38  .1  =  1  ,N 
.I1=J1+J 
AMI J)=A( J1 I 
38  CONTINUE 

WRITEI3,34HAM(.I).J=1,N) 

34  FORMAT!  • 36  I  GENVALUE S ' / ( 1 X , 5E 14 .7 ) ) 

41  WRITE(2,3ft)(AM(  J),.I=1,N) 

3ft  FnRMAT(5E14.7) 

NM=N*N 

WRirE(3,3?)|(P(  I  ).l  =  l,NN) 

35  FORMAT! ' OE I GENVECTUR 5 ' / ! I X .5E 14 .7 1 > 

42  WR I TEI2. 3ft) IP!  I)  .  I  =  l.NN) 

N1=N*!N+1 )/2 
WRITE!3,2ft) !A! J),J=1,N1  I 

2ft  FORMAT! 'ODI AGONAL  MATR I  X ' / I  1 X 1 7E1 1  .4 )) 

.11=0 

DO  28  J  =  1  ,N 
AM! J)=0. 
no  29  I  =1  ,N 
111  =0. 

.12=  I 

00  33  K=1,N 
.I3  =  J1-H< 

U1=U1+TT(.I2)*P!.I3) 

J2  =  J2-fN 
33  CONTINUE 

.|4  =  J14.I 

Sl=Ul*r.i:N.lG!P!,)4)  ) 

AM!  J  )  =  AM  J  )  +  Sl 
29  CONTINUE 
J1=J1+N 
28  CONTINUE 

WRI  TE  !  3. 38  )  !  AM!  J  I  ,.I=1,M) 

38  FORMAT!  'ORAYLEIGH  OUOT I  ENT S • / ( IX t5E 14.7  I ) 

STOP 

END 

JOATA 

10  3ft  0.2808^826+00 

0.0000  5.3f08  10.ft048  15.5027  18.?ft48  18.2548  15.5027  10.ft048  5.3ftOR  O.OOOf 
-10.0000  -8.5-20  -8.4705  -5.2884  -2.5558  2.5558  ft. 2884  8.4705  8.5320  10.0000 
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5.0000  15.0000  25.0000  35.0000  55.0000  55.0000  55.0000  75.0000  S5.0000  ^5.0000 
105.0000115.0000125.00001  35 . 0000155 .00001 55 .0000 1 5  5. 00001 75. 0000 1  H5. 0000 1 95 .0000 
205. 00002 15. 0000225. 000023 5. 0000255. 0000255. 0000255. 0000275. 00002H5. 0000295. 0000 
305.0000315.0000325.0000335.0000355.0000355.0000 
iSTOP 
/  ♦ 

// 

PRINTED  OUTPUT 

N  M  HK 

10  35  0.290flP82E+00 


X 

0.0000  5.3508  10.5058  15.5027  19.2558  19.2558  15.5027  10.5058  5.3508  0.0000 

Y 

-10.0000  -9.5320  -8.5705  -5.2895  -2.5559  2.5559  5.2895  8.5705  9.5320  10.0000 

PHI 

5.0000  15.0000  25.0000  35.0000  55.0000  55.0000  55.0000  75.0000  85.0000  95.0000 
105.0000115.0000125.0000135.0000155.0000155.0000155.0000175.0000185.0000195.0000 
205.0000215.0000225.0000235.0000255.0000255.0000255.0000275.0000285.0000295.0000 
305.0000315.0000325.0000335.0000355.0000355.0000 


EIOENViLUES 

0.9557155F  +  02  0.8078329F  +  02  0.519531  3F  +  02  0. 555859  7F  t-02  0.2875992E+02 
0.2551213E  +  02  0.9539283F-t-01  0. 32755  1  2F  +  0 1  0.55075571-00  0.  1  528523F+00 

FIGENVFCTOHS 

0.2903522E  00  0 .2 30562 1 E-05  0.5855855E-01  0 .525 1 972F -08-0. 355752 1 E  00 
-0.9729786F-07-0. 588551  IF  00-0.  2973702F-07-0. 19B20P3F  00  0.5377555E-05 
0.19R2075E  00  0.2557125F-06  0.58H5510E  00-0. 309 1 38 7F -05  0.3557533E  00 
-0.1850908E-05-0.6855RR5F-01  0.  27 1 355  7E-05-0. 2903550  8  00  0. 5529575E-05 
-0.89159558-01  0.2151319F-05  0.2720950F  00-0.15330208-05  0.33087878  00 
-0.99515518-07-0.1  1726598  00  0. 38587 1 5F-P5-0. 55  20 352F  00  0.29575588-05 
-0.55293788  00-0.53687208-07-0.11725738  00-0.5255318^-05  0.33087B0E  00 
0.30878538-05  0.27209518  00  0.  7309595E-05-0.  89 1 595  3  8-0 1  0.21002978-05 
0.529R515E  00  0.527025HE-05  0.530R088E  00  0.15739078-05  0.29205238-01 
0.95195298-05-0.77851998-01-0.55255188-05  0.15351528  00-0.88555358-05 
0.15351508  00  0.17715758-05-0.77852058-01  0.55551558-05  0.29213558-01 
-0.6527H558-05  0.53082038  00-0.  15295588-05  0.52985558  00-0. 752 1  >5978-05 
-0.10A7320E  00-0. 19358718-05-0. 5557035E  00-0.1  159  3898-05-0.395  85588  00 
0.1512510E-06  0.88930598-01  0.55750798-05  0.15555358  00-0. 72591 588-05 
-0.16555556  00-0 .8 1 52 255F-05-0  .HH9?91 8F-01  0.52928558-05  0.39585508  00 
0.6803169E-06  0.55569578  00-0. 58201978-06  0.10572218  00-0. 1 509 1 32E-05 
-0. 377 35  77 E-0 1-0. 2 795  0 798 -05  0 , 1 8755058 -01 -0 .50975098  -  05-0 .5 157685E  00 
0.35259278-05-0.55731318  00  0.23555558-05-0.11802508  00-0.15773338-05 
-0.11801706  no  0.32595738-05-0.55730588  00  0.55751258  -05-0.51677228  00 
0.22525838-05  0.18758988-01  0.25139078-05-0.37723288-01  0,32589598-05 
-0.53597128  00-0.35588558-05-0.1  1951078  00-0 . 8S0585 7(  -05  0.10509198  00 
-0. 191 82 2 18-05-0. 28  1  OOH 18  00-0.5109  35  38-05-0.  32897 38 8  00-0.271 1 8Q98-05 
0.32R9759E  00  0.19851958-05  0.28102018  00-0 . 1 35 7 1 35 r -05-0 . 1 0508 1 18  00 
-0.32871228-05  0.11951028  00  0.55172788-05  0.53595928  00  0.31988258-05 
-0,51397688  00-0.92519598-06  0.W210958  00  0.35535568-06  0.20250038  00 
0.19255788-06-0.19551058  00-0.97627018-06  0.  389957  8  00  0. 698 7 1  3 38-06 
0.3899392E  00-0.39036928-05-0.  9690798  00  0.529599  8-07  0.20290228  00 
0.20553038-05  0.31210528  00  0.19039588-05-0.91397598  00-0.25853288-05 
0.32109838  00  0.91926978-05-0.39951928  00-0.95519538-05  0.23552598  00 
0.225''7978-05  0.19750038  00-0.99599518-05-0,99753998  00-0.91580958-07 
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0.«^7f.404E  00-0.4‘»<44?04E-0(S-0.147^034E  00-0 . 3?Sflft‘j4F -04-0.  ?364?09E  00 
0.2H7HH05e-06  0.3446209F  00-0. 54<»633?F-04-0. 3211 0?qE  00  0. 5P 1 7425E-04 
0.1<Jft22‘54E  00  0.1733023F-05-0.3774212E  00-0 . 3023024F -05  0.41fl3334E  00 
0.37l0524E-05-0.3425n54F  00-0.  3 12 1 30<)F-05  0.m7757F  00  0.  1 42  5577F-06 
0.1117430E  00-0.427^':  1  4F-04-0.3424925F  00-0 . 027MR49F -07  0.41H322(SF  00 
0.l33fl545E-04-0.377hl49E  00-0. 35«4534F -04  0.1942231F  00  0. 1 24 1 534F-07 
-0.119H015F  00-0.1h55H49F-04  0.2539174F  00  0 .4 1 1 1 4« 7F -04-0 . 379 1949E  00 
-0.35I4473F-04  0.3907340E  00-0. 3590495E-04-0. 353 U 59F  00  0. H053 73 1 F-OF 
0.3531210F  00-0.23H4469F-05-0.39074fl4F  00  0. 294H/.22F-05  0.37920fl4e  00 
-().35245nF-05-0.25392/9F  00  0. 2444993F-05  0.n9M04HF  00-0.  U44047F-05 

n I  AGONAL  M/T«IX 

0.9447F  +  0;'  n.OOOOE  +  OO  0.fi542F-08  0.32836-10  0.H078E  +  02  0. OOOOE  +  00-0. 57R 1 6 -07 
-0.2912F-1'  0.21886-05  0.3970E-10  0.4195F+02  O.OOOOE+00  0.50326-05  0.10R4F-08 
-0.27376-0M  0. 55816-11-0. 2503E-04  0.10736-10  0.544RF+02  O.OOOOE  +  00-0.  173HE-07 
0.L417F-1-  0. 32346-06-0. 4128E-10-0. 82526-05  0. 509  7F-09-0 . 8921 6-04  0.4R99F-09 
0.28746+0;'  0.00006+00-0.81406-05-0.54996-08-0.15046-05  0.1549E-09  0.2023E-04 
-0.21996-1-0.35326-07  0 . 6/966-1 1-0. 1455F- 1 0  0.00006+00  0. 25516+02  0.00006+00 
0.14556-10-0.24876-13  0.31496-07-0.31716-11  0.15446-09  0.12876-13  0.4015F-09 
0.12046-1  0.53946-05-0.45456-09-0.75476-09-0.43016-14  0.9439F+01  O.OOOOF+OG 

-0.20626-09-0. 4 1006-1 2-0. 291 36-10  0.22296-15-0.35786-07-0.13776-11-0.  13406-10 
-0.54936-1;.  0.28906-05-0.88686-10  0.25496-05-0.14506-09  0.00006+00-0.11106-14 
0.32756+0.  0.00006+00  0.37136-05  0.24746-09-0 . 25046-05-0. 88746 -10-C . 7 1 84F -0 7 
-0.14776-10-0.23716-09-0.17016-12  0.34206-09  0.40806-13-0.44146-04  0.33446-0'' 
0.53906-09-0.65326-12-0.27826-07  0.43716-10  0.55086+00  0. 00006+00-0 . 288 IF-O-) 
0.23506-1  -0.11  196-07  0.  1 2 1 3F-1 1-0. 1 3476- 10  0. 338 36-1 4-0. 4 5776-05  0.2410F-nc 
-0.23286-09  0.1  5546-14  0.79326-09  0.53116-14-0.45336-05  0.7744E-09-0.12346-0‘- 
-0.34536-i;'-0. 14556-10-0.64416-15  0.14286  +  00  0.00006+00 

HAYLFIGH  OHflTIENTS 

0.9447105  :  +  n2  0.80783076  +  02  0.41953146  +  02  0.54489946+02  0.28  740136+02 
0.2551244i:+02  0.96393046+01  0.32754006+01  0.55074946+000.14283916+00 
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V.  PATTERN  SYNTHESIS  WITH  CONSTRAINED  SOURCE  NORM 


The  program  of  this  section  Implements  the  theory  of  section  V  of  Re¬ 
port  Ho.  2.  Punched  card  Input  data  is  read  Into  the  program  according  to 

READ(1,30)  N.M,  NO,  BK,  EPS 
30  FORMAT  (313,  2EU.7) 

READd,  10)  (X{I>,  l-l .  N) 

10  FORMAT (10F8. 4) 

RFj\n(l,10)  (Y(l),  I-l,  N) 

READd, 10)  (PHI  (1),  1«1,  M) 

READd, 23)  (AM(J),  J-1,  N) 

23  FORMAT (5E14, 7) 

NN  -  N*N 

READd, 23)  (P(J),  J-1,  NN) 

DO  22  JG-1,  NG 
READd, 23)  CN 
READd, 10)  (G(I),  1=1,  M) 

23  CONTINUE 

The  X  and  y  coordinates  of  N  point  sources  are  read  In  as  X  and  Y.  The 
specified  pattern  g  of  (2-16)  Is  read  in  as  G  at  M  angles  PHI,  The 

^  'N/ 

eigenvalues  AM  and  eigenvectors  P  of  the  matrix  lT*T]  of  (2-35)  have  been 
punched  out  by  the  program  of  section  IV  of  this  report.  The  constraint  C 
of  (2-24)  is  read  in  as  CN.  The  iterative  process  for  finding  the  root  of 
(2-45)  is  terminated  as  soon  as  F(a)  _  EPS*CN.  DO  loop  22  makes  it  possible 
to  deal  successively  with  NG  different  specified  patterns  at  the  same  propa¬ 
gation  constant  BK. 

Minimum  allocations  are  given  by 

COMPLEX  P(N*N),  T(N*M),  G(M),  TG(N),  C(N), 

AI.P(N),  F(N),  TF(M) 

DIMENSION  X(N),  Y(N),  PHI(M),  AM(N),  AM2(N), 

CC(N),  TFM(M),  FP(N),  FM(N) 
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Nested  DO  loops  14  and  15  store  the  matrix  element  given  by  (2-18) 
in  T((m-l)*N+n) .  IX)  loop  25  stores  in  and  replaces  G  by  its 

complex  conjugate.  DO  loop  20  stores  the  complex  conjugate  of  [T*]g^  of 
(2-4.!)  In  TG.  Nested  DO  loops  26  and  27  store  Cj  of  (2-42)  In  C(J)  and 
!c,|  ’  in  CC(J). 

DO  loop  31  stores  F(0)  of  (2-45)  In  S6.  If  S6  •  0,  DO  loop  68  carries 
out  Newton's  method  starting  with  a  equal  to  the  center  eigenvalue.  Refer¬ 
ring  to  expression  (2-45),  a  is  stored  in  both  AL  and  AS  while  F(a)  Is 
storod  in  FS.  Statement  60  calculates 

F(a  ) 

“j+1  "  “  “j 

If  a  ,  <0  then  a...  is  set  equal  to  zero.  The  subscripted  a's  appearing 
'+1  “  J+l 

in  (  .)  are  not  to  be  confused  with  those  of  (2-40).  DO  loop  91  puts  of 
(2-40)  in  ALP(l). 

DO  loop  33  stores  f  of  (2-41)  in  F  and  accumulates  i|?lP  in  FF.  Upon 
exit  from  DO  loop  33,  F(J2)  will  be  the  element  of  1  largest  in  magnitude. 
DO  1 )Op  92  stores  the  synthesized  pattern  T?  of  (2-23)  in  TF  and  the  magni¬ 
tude  i  of  the  elements  [T]f  in  TFM.  ||f|p  is  accumulated  in  TFS  and  c  of 
(2-2J)  in  E.  Just  after  DO  loop  92,  Q  defined  by  (2-15)  is  stored  in  Q  and 
E  is  normalized  by  dividing  by  GG.  Statement  59  punches  f  on  cards  for  pos- 
slbl'!  input  into  the  program  of  section  IX  of  this  report.  DO  loop  41 
normillzes  the  element  of  f  largest  in  magnitude  to  unity  and  then  stores 
the  magnitudes  and  phases  in  degrees  of  f  in  FM  and  FP  respectively. 

According  to  the  printed  outputs  of  the  pattern  synthesis  with  con¬ 
strained  source  norm  program  of  this  section  and  the  pattern  synthesis 
program  of  section  II  of  this  report,  it  has  been  possible  to  reduce  the 
source  norm  squared  from  13.366  to  4.0  with  only  a  very  slight  increase 
in  t  le  pattern  synthesis  error. 


F(a,) 


1^1 


(1) 
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LISTING  OF  PATTERN  SYNTHESIS  WITH  CONSTRAINED  SOURCE  NORM  PROGRAM 


//  (n034,FF.  lSS.li4l,'MAIiT7iJOE*  iREr,ION»200K 

//  EXEC  WATEIV 

//G0.ET02F001  1)1)  S Y SOUT  =  H . OCR*  ( RECFMiE  ,RLK S I  7 E  =  RO  > 

//GO. SYS  IN  1)1)  * 

$JOB  MAUT7  .T  IME=1  .PAGES=40 

COMPLEX  U,(J1,C(JNJG,P(  I  001 .  T  (  3/.0 )  ,0  (  3ft  )  .  TG  (  10  ) .  C  (  10 1 .  AL  P  (  10  I 
COMPLEX  1)2, F  (  101  .TF(  3ft» 

01  MENS  I  ON  X  (10)  ,Y  (  10)  ,PHl  (3ft»,AMnO)  .AM2(  10)  ,CC(  10l,AS(3n 
DIMENSION  FS (30)  ,TFM  (3ft)  .FP( 10) .EMI  10) 

P!=3.1ftl5<)3 
CP=P1 /IflO. 

0=(0..1.) 

RFA0(1,30)  N.M.NG.HK.FPS 
30  FORMAT! 313.?Elft.7) 

WR1TEI3,32)  N.M.NG.HK.FPS 

32  FORMAK'O  N  M  NG  •  ,  ftX  , ‘BK*  .  1 1  X  .  •  E  PS  •  / 1  X  ,  3 1  3. 2E  14. 7  ) 

REAOl  1 ,10)  (X  (  I  )  ,  I  ==1  ,N) 

10  FOKMATdOFH.ft) 

WRITE(3.13)(X( I  )  ,1=1  ,N) 

13  FORMAT! 'OX '/ ( 1 X , 10FH.4) ) 

RFA1)(  1 ,10)  (Y  (  I  )  ,  I  =1  ,N) 

WRITE(3,3/)(Y( I  ),I=1,N) 

37  FORMATC 'OY •/ ( 1 X, lOFfl.ft) ) 

READ! 1 , 10) (PHI  ( I  ) , 1 =1  ,M) 

WRI  TF(3,1  1 )  (  PHI ( ( ) , I =1  ,M) 

11  FORMAT! 'OPHI •/(lX,10Ffl, 4)) 

REAOl  1  ,23)  ( AM(  .1  )  ,>l  =  l  ,N) 

23  F0RMAT(‘jE)4.7  I 
WRITE(3,24|(AM( J),J=1,N) 

24  FORMAT! 'OAM</( IX, SF14.7)) 

NNaN*N 

7M»M 

HEA0(1,23)(P(.II,J=1,NN) 

WR  I  TE  (  3,21  I  (P  (  J  )  ,.l  =  l  ,NN) 

21  format  ( '01"/ (  IX  ,  SF  14.7)  ) 

no  24  .1=1, N 

AM2( J)=AM( J)mAM( J) 

24  CONTINUE 
J1=0 

no  14  K=1,M 
S1=PHI  (K)«=CP 

cs=cas( SI ) 

SN=S1N(S1 ) 

DO  15  1=1, N 
J1=J1+1 

S1=RK*(X(I  )4r.S+Y(  I  )*SN) 

T(ji  )=cns(si )-ui*siN(si) 

15  continue 

14  continue 

on  22  .IG  =  1,NG 
RFAn(1.23)  CN 
WRITE(3,2R)  CN 
28  FORMAT! 'OCN=' ,E14.7) 

READ! 1 , 10) (G( I  )  ,  I  =  1 ,M  ) 

WRI TE( 3,12) (G( I  )  ,  I =1  ,M) 

12  FORMAT ( 'OG '/( 1 X , 10EH.4 ) ) 

GG  =  0. 

DO  25  .)=1,M 
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iii3CiJN.i('.(r>(  j )  I 
r,G=Gr,+G(  j)»ui 
G(J(=U1 

25  CGNTINUe 

011  14  J  =  l,N 
J2  =  J 

TGI J (=0. 
on  20  K=1,M 
TGIJlrTGIJl+TI J?»*GIK( 

.I2  =  J2+N 
20  CONTINIIR 
14  CflNTlMlIF 
.11=0 

on  25  .1=1, N 

01=0. 

Dll  27  1  =  1, N 
J1=J1+1 

111  =U1+P(J1  )  fTGI  I  ) 

27  CONTINOt 

C I J l=CnNjG( U1  ) 

CC(J)=t)l*C(.l) 

26  CONTINUE 

W«ITe(3,200)  (Cr,(.l),.l  =  l,N) 

200  KIHMATI  <OCC'/(  IX.bFH.TI) 

.IAM=  (N+1  1/2 

f  P  =  ePS*CN 

201  S6=CN 

on  31  K=l,N 
S6  =  S6-r.r.(K  l/*M2tK  I 
31  CONTINUH 
.14=1 

IF(S6)  66,65,65 

65  *L=0. 

AS(J4)=4L 
FS(J4|=S6 
Gil  TO  H7 

66  AL=AM(JAM| 

AS ( J4  )=AL 

on  6H  ,1  =  1  ,30 

S2=CN 

S5  =  0. 

no  67  K=1,N 

S3=AM(K  )+4L 

SA=S3*S3 

S2  =  S2-Cr.(K)/S6 

S5  =  S5-fCC  IK  I  /  (  S4*S3) 

67  CONTINUF 
FS  I  .14  l  =  S7 

IF(AFS(F.SI.I9)  l.LF.FP)  GO  TO  87 
6f  AL=AS(.|4)-.5't‘FS(J4|/S5 
IFIAL.LF.0.1  Al  =0. 

J4  =  .|4+l 
AS! J4)=AL 
6f  CONTINUE 

WKI  TK  3,214  ) 

21  ‘  FORMAT  I ' ONFWTONS  METHOD  DIO  NOT  CONVEBOF') 
8’  no  41  K  =  1  ,N 

ALPIK )=C IK ) / ( AM  IK  l  +  AL 1 
41  CIINTIMIIF 

WW  I  TE  (  3,40  I  I  AS  I  K  )  .K  =  l  ,.I4) 
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no  FIIRMATC  convergence  f)E  LAGRANGE  MOLT  I  PL  t  ER  •  /  (  I X 1 7E 1 1  .A )  ) 
WRITE(3,100)(ES(KI,K*l,jn» 

100  FflRMATC  EUNCTION'/(  1X,7EU.4)  » 

FF  =  0. 

S2  =  0. 

Of)  33  I  =  1,N 
F(I )=0. 

J4=I 

DO  44  Krl ,N 

F(  I  )  =  F(I  l■*•ALP(K)♦P(  J4| 

J4=J4+N 

44  CONTINUE 

51  =  F ( I (♦CONJOIF (  I  I  I 
FF=FF+S1 

IFISl.LT.S?)  GIJ  TO  33 

52  =  S1 

J2=I 

33  CONTINUE 
E=0. 

TFS=0. 

J1=0 

on  92  KK  =  1  ,M 

U1=0. 

on  93  K=1,N 
J1=J1+1 

U1«U1+T( J1 )*F(K) 

93  CONTINUE 
TF (KK )=U1 
U2=C0NJG(U1 ) 

Sl=Ul*U? 

TFS=TFS+Sl 
TFMIKKIsSORTISn 
U1=U2-G(KK ) 

S1=U14C0NJG(U1 ) 

E=E+S1 
92  CONTINUE 
q=ZM*FF/TFS 
E=E/GG 

WRITE(3,48)  GG.FF.EtO 

48  FORMAT! '0GG=' ,F14.7,'  FF=',E14.7,*  E»'.E14,7,»  0*',E14.7) 
WRITE(3,4^))(F(  I  )  .1=1  ,N) 

45  FORMAT! '0F'/( IX, 5E14. 71) 

59  WRITE (2,23) (F(I ) , I =1  ,NI 

Ul=l./F(.)2) 

00  41  J=1 ,N 
F(  J)=U14F(.J  I 
FM(J )=CARS(F ( J ) 1 
S1=REAL(F(J) ) 

S2  =  AIMAr,!F(  J  )  ) 

Fi'lJ)  =  (ATAN2(S2.Sl  )  )/CP 

41  CONTINUE 

WRITt(3,42)  (FM(J  »  ,.l  =  l  ,N) 

42  FnRMAT('0FM'/(lX,10F7.3M 
WRITF(3,43) (FP( J),J=1,N) 

43  F0RMATI'0FP'/(1X,10F7.1I) 

WRITE(3,49)(TF(K),K=1,M) 

49  FORMAT! 'OTF'/dX.lOFH. 41) 

WRITb!3,55) !TFM!K  )  ,K  =  1,M) 

56  FORMAT! 'OTFM'/! 1X,10FH.4)) 

22  CONTINUE 
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0.3?483|3r  00  0.4770?4ft8-03  0.44080888  00  0.14734078-03  0 . ?4?03?38-0 1 
0 .431  43?4t--03  -O.  7  7  8 4  1  448 -0  !  - 0 . 3  3  ?4 3  1  ft  -  03  0.1  >331  3?8  00-0.  ft  84  33 338 -03 
0.13331308  00  n.  I  7  71  3 /3I--03-0.  7 /H4?0'.3-01  0.43  34)348-03  0  .  ?4?  1  3338 -0  1 
-0.7.3?  7X338-03  0.4  0X20-0  Oil-0 .  1  4  ?4  34  8  8 -03  0.3?48343F  00-0 . 7  3  ?  1  34  78 -03 
-0.104  7  3208  on  -0. 1  4338  7 1  8-03-0.34  3  70338  00-0 .  1 1  33)  ))fl4r -03 -0 . 344  8  33  ft8  00 
0.141?3108-03  0.88430448-01  0 . 44 7 30 74f -03  0.13333338  00-0 ,  7 ?4 4 1  3»F -03 
-0.13333331-  no -0.8  1  3??34(  -(i/,-0.884.'>4|(  8-0!  0 . 3?4 ?« 34 8 -03  0.34483308  00 
0.38031  ()48-n7-,  0.343343  7>-  00-0.3  8  2014  78-03  n.l047??lF  00-0.  1  304  1  328-03 
-0 ,37  7  34778-01  -0 .2  /440  .’48-  43  0.1  M  7  3  038-0! -0 .404  73048-03-0.41373848  00 
0.332  34?  ,'r-n3--n.  33  7  3  !  -il  >  00  0 .?  3  •>  33348 -03-0 .  1  1  x  02308  00-0.  14  7  733  38-03 

-0,11  801  708  on  0.32344  7'<[ -03-0,33730388  00  (' .  3  3  74  |  ?48 -03 -0 .4  1 3  7  7  ??8  00 
0.224248 28-03  0. lx/xx4XX-n!  0 . ?4  1  3407 8-03-0.  3 7 7 ?3 ?8 8-0 1  0 .  3?3 ft 44 48 -03 

-0.3334  7  128  on  -0.33488338-03-0. 1  143107F  00-0 . 8 304  84 78 -03  0.1  0304148  00 
-0.141 82218-03- 0,28  I  008  18  00-0. 3 1 04  34  38-03-0. 328  4  7  3X8  00-0.271  1  8  448-03 
0.32847348  00  0,1483)4x8-03  0.27(102018  00-0.13371338-03-0.10308118  00 
-0 , 328  7  I  228-03  0,11431028  HO  0.  33  1  /2788-03  0. ’■3343428  00  0.  3  ]  4RH238-03 
-0.41347388  nn-n. 423 14448-03  0.31210438  00  0.33333438-03  0.20240038  00 
0 . 14?  S4  7718-03- 0. 1434  I  038  00-0.4  7  3  27018-03  0.  3R4P4V78  00  0.  34  8  71  338-03 
0.38443428  00-0,34033428-03-0. 14340748  00  0.32444418-07  0,?0?40??8  OQ 
0.20333038-07-  n.31?ir3?8  00  0.  140343X8-03-0.41347548  00- 0. 23 ft 33288-03 
0.  ’-2  1  04X38  no  0, 414234  78-03-0. ft44314?8  00-0.4  33  143  38-03  0.23332348  00 
0.2237 747x-nx  0.1473nn’7  00-0.44440318-03-0.44733448  00-0.41380438-07 
0.44734048  00--n. 44442048-03-0. 1-730148  00-0 . 32  X  X334  8 -03 -0 . 2  3  33  2048  00 
0.2X78X038-03  0.34432048  00-0 . '> 3433 328 -03-0.  1 2 1  1 0 ?48  00  0.34  173238-03 
0. 14322348  no  0.17 170238-03-0.2  7  -32128  00-0.30230248-03  0.417  33348  00 
O.371O3248-n3--0.332‘'nS48  00-0.31  2  13048-03  0.  11177378  00  0.  14235778-03 
0.11173208  nn-0. 42/3417, 8-03-0. 2324423I-  00-0  .  X  ?  7M  8  341  -07  0.41832238  00 
0 . 1  338343*  -07>- 0. 37  7  31  448  00-0 . 2 <- x4 3 33(- -03  0.  1  43 ? ? 3 1  8  00  0.  1  24  1  3348-07 
-0.1  1*8X0131  nn-n.  1333X448-03  O. 23341748  00  0 .4  1  1  |  4R  78 -03 -0 . 3  74  1 4448  00 
-0,231  34  738-03  0.34073308  00-0 . 3  340443  8 -03-0 .  33  3  1  1  34(-  00  0.80337318-03 
0.33312101  00--n. 23844348-05-0.34074838  00  0 . 7>J3X4  228  -  03  0.37420R38  00 
-0 .332 33  1  18-03-  0. 23  342  74P  0  0  0.2  4  4  4  4  4  3  8  -  0  3  0.  1  1  4  X  0  3  8  8  0  0-  0.  1  14  3  0  4  78  -  03 
0 .40000008-701 


3.0141 

o.noon 

2.8317 

o.noon 

2 . 2337 

n.ooi'n 

1  .  74  24 

n.onoo 

1.414? 

0.0000 

1  .??()H 

o.onno 

1 . 1 024 

o.rooo 

1 .0331 

o.nono 

1 .on2p 

0.0000 

0.0000 

o.nono 

0.0000 

o.noon 

o.onno 

n.nooo 

n.ooon 

n.ooon 

0.0000 

n.ooon 

0.0000 

o.nono 

n.onoo 

0.0000 

0 . noon 

0,0000 

n.oono 

n.onoo 

o.noon 

o.onno 

0,0000 

n.onoo 

0 . on on 

n.nooo 

0 . noon 

o.onno 

o.onno 

o.onno 

0. 0000 

0. 0000 

o.nono 

n.ooon 

0 .0000 

o.noon 

0.0000 

n. oonn 

o.'inf.n 

o.onno 

0.0000 

n.oono 

o.noon 

n.nooo 

o.noor 

o.noon 

o.onno 

0 .  n  '<0  0 

().0f,00 
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0.0000 

0, noon 

0.0000 

0. oonn 

n.onoo 

0.0000  0.0000 
ISTOP 
/« 

// 

PRINri-O  OUTPUT 

N  M  NO  RK  PPS 

10  3A  1  0.?')0HHR;(-  +  00  O.lOOOOOOF-02 


X 

0.0000  S.36PR  10.ft04H  lP.?ft4R  19.2iS4R  15.3027  10.(S04R  5.360R  0.0000 

Y 

-10.0000  -9.5320  -H.4705  -5.2H94  -2.5559  2.5559  5.2R94  H.4705  9.5320  10.0000 

PHI 

5.0000  15.0000  25.0000  35.0000  45.0000  55.0000  55.0000  75.0000  fl5.0000  95.0000 
10  5.00001  15.00001 25.0000135.000014  5.0000155.0000155.00001  75. 0000 1 fl5 . 0000195 .0000 
205. 000  02 15. 0000225. 000023 5. 0000245. 0000255. 0000255. 0000275. 00002fl5. 0000295. 0000 
305.0000315 .0000325.00003  35.0000345.0000355.0000 

AM 

0.9457|*.51-  02  0.HO7H329F  02  0.5195313F  02  0.545R497F  02  0.2R75992F  02 
0.2551213F  02  0.9n342H3F  01  0.3275412F  01  0.5507557F  00  0.152R423E  00 

P 

0.29035226+00  0 .2 30452 1 F-05  0.5R44R55F-01  0 .424 1 9 726 -OR-0 . 354 74 2 1 F +00 

PLUS  39  muhf  linfs  uf  pkintfo  output 

CN=  0.4000000F+01 
0 


5.0191 

0.0000 

3.8537 

0.0000 

2.3552 

0.0000 

1.7434 

0.0000 

1 .4142 

0.0000 

1 .2208 

0.0000 

1.1034 

0.0000 

1.0353 

0.0000 

1.0038 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.oono 

0.0000 

O.OOOO 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

CC 

0.25525516  +  03  0.42031546  +  03  0 . 2 83405 1 F  +  02  0, 2705 12R6 +03  0.35131156+03 
0.131  19166+03  0.1149905F+02  0.  12501 256  +  02  0. 2570739F  +  00  0.27551406  +  00 
CONVFPr.FNCF  fIF  LAGRANOF  MULTPLIEH 

0.28756+02  O.OOOOF+00  0.71236-01  0.15346+00  0.21906+00  0.24406+00  0.24546+00 
FUMCT ION 

0.3  7316  +  01-0.93556  +  01-0.37  556  +  01-0. 1  287F  +  01 -0. 29 1  5F+00-0. 2 370F-01 -0 . 1 84 lF-03 
r,G=  0.55545096  +  02  6F=  0.40001 40F  +  0 1  F=  0.324  34236  +  00  0=  0.404P585F  +  01 


6 

-0.34271296  00  0.23984996  00  0.79558796  00  0.1R85415F  00-0.93774596  00 
-0.33021876  00  0.15213596  00  0.78945996  00-0. 409552 7 F-03-0.  2 1 25 1  886  00 
-0.30357586  00  0.4592439F  00-0 . 1 559025F-01-0. 1 82 1  1976  00-0.19531356  00 
0.51733516  00-0.5551  524F  00-0.43237236  00  0.4281348F  00-0.  1 7 1  7509F  00 

6M 

0.421  0.822  1.000  0.811  0.214  0.552  0.184  0.552  0.708  0.454 
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-54.4 

173.9 

0.0  -121.0  70 

.5  -74.5  45.4 

-91.8 

18.5 

138.7 

TF 

2  ,4114 

0.0952 

2.9077 

-0.0472 

2.4434  -0.3090 

1.7975 

-0.3829 

0.9121 

-0.0844 

0.4><27 

0.3438 

0.5554 

0.5280 

0.4874  0.2744 

0.6985 

-0.0894 

0.4039 

-0.2834 

0.2M41 

-0.3522 

0.2310 

-0.4421 

0.0845  -0.5259 

-0.1097 

-0.5108 

-0.2429 

-0.4192 

-0.25<*‘» 

-0.3459 

-0.1844 

-0.3355 

-0.0742  -0.3590 

-0.0032 

-0.3730 

-0.0323 

-0.3757 

-0.1494 

-0,4062 

-0.2566 

-0.4842 

-0.2449  -0.5705 

-0.1000 

-0.5771 

0.1083 

-0.4517 

0.2513 

-0.2044 

0.2142 

0.0454 

0.0097  0.2043 

-0.2235 

0. 1424 

-0.3873 

0.0718 

-0.5540 

0.0555 

-0.7741 

0.0575 

-0.8721  -0.0190 

-0.5454 

-0.1191 

0.2513 

-0.1057 

1 .3812 

0.0;>43 

TFM 

2  .4135 

2.9081 

2.6414 

1  .8379 

0.9140  0.4125 

0.7445 

0.7411 

0.4052 

0.4934 

0.4401 

0.49flft 

0.5329 

0.5225 

0.4845  0.4324 

0.3828 

0.3444 

0.3731 

0.3771 

0.4329 

0.5't97 

0.4217 

0.5857 

0.4445  0.3240 

0.2259 

0.2044 

0.2743 

0.3939 

0.5547 

0.7'42 

0.8723 

0.57HO 

0.2727  1.3814 
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VI.  MAGNITUDE  PATTERN  SYNTHESIS  WITH  CONSTRAINED  NORM 


The  program  of  this  section  Implements  the  theory  of  scctlone  IV  and 
V,  Report  No.  2. 

Punched  card  data  Is  read  In  according  to 

READ(1,30)  N,M,NG,N9,BK,EPS 
30  F0RMAT(4I3,  2E14.7) 

READ(1,10)(X(I),  I-l,  N) 

10  FORMAT (10F8. 4) 

READ(1,10)(Y(I),  I-l,  N) 

READ (1,10) (PHI (I),  I-r,  M) 

READ(1,23)(AM(J),  J-1,  N) 

23  FORMAT(5E14.7) 

NN  -  N*N 

READ(1,23)(P(J),  J-1,  NN) 

DO  22  JG-1,  NG 
READ(1,23)  CN 
READ(1,10)(G(I),  I-l,  M) 

22  CONTINUE 

The  X  end  y  coordinates  of  N  point  sources  are  read  In  as  X  and  Y,  The 
eigenvalues  AN  and  eigenvectors  F  of  the  matrix  [T*T]  of  (2-35)  have  been 
punched  out  by  the  program  of  section  IV  of  this  report.  The  starting 

value  of  h  e  In  (2-21)  Is  read  In  as  G  at  M  angles  PHI  In  degrees  In 

B 

the  plane  of  the  point  sources.  The  magnitude  of  G(m)  Is  the  magnitude  of 

j  g 

the  specified  pattern  at  the  m*^^  angle.  The  phase  portion  e  of  G(m) 
Implicitly  gives  the  starting  value  for  DO  loop  22  makes  It  possible 
to  deal  successively  with  NG  different  specified  pattema  at  the  same  propa¬ 
gation  constant  BK.  The  constraint  C  of  (2-24)  Is  read  In  as  CN.  The 
Iterative  process  for  finding  the  root  of  (2-45)  Is  terminated  as  soon  as 
F(a)  <_  EPS*CN.  Step  2  of  the  Iteration  procedure  described  In  section  IV, 
Report  No.  2,  la  repeatedly  executed  until  the  pattern  synthesis  error  e  of 
(2-19)  ceases  to  decrease,  but  not  more  than  a  maximum  of  N9  times. 
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f 


Hlnlmum  allocations  are  given  by 

COMPLEX  P(N*N),  T(N),  TP(N*M),  G(M),  C(N), 

ALP(N),  TF(M).  F(N) 

DIMENSION  X(N),  Y(N),  PHI(M),  AM(N), 

AM2(N),  H(M),  FF(N9+1),  E(N9+1), 

Q(N9+1),  CC(N),  TFM(M),  FM(N) ,  FP(N) 

X)  loop  15  stores  the  row  of  the  matrix  [T]  of(2-18)  In  T.  DO 
loop  I'j  stores  the  K^^element  of  [T^j  In  TP( (J-1)*M+K) .  DO  loop  25  stores 
h  of  (2-19)  In  H,  accumulates  Ilg^jlF  of  (2-14)  in  GG  and  replaces  G  by  Its 
complet  conjugate. 

Steps  1,  2,  and  3  of  the  iteration  procedure  described  in  Section  IV, 
Report  No,  2,  are  executed  for  the  K9  time  inside  DO  loop  19.  Dl'  loop  18 
stores  Cj  of  (2-42)  In  C(J)  and  ICjp  in  CC(J).  DO  loop  31  stores  F(0)  of 
(2-45)  In  S6.  If  S6  <  0  DO  loop  68  carries  out  Newton's  method  starting 
with  o  equal  to  the  root  of  (2-45)  for  the  previous  value  of  K9.  a  Is  stored 
in  both  AL  and  AS(J9)  while  F(a)  of  (2-45)  resides  in  FS(J9).  Statement  60 
calculates  (1), 

If  £  0  then  equal  to  rero.  The  subscripted  a's 

appearing  in  Equation  (1)  are  not  to  be'  confused  with  those  of  (2-40). 

DO  loop  91  puts  of  (2-40)  in  ALP(i)  and  accumulates  ||f  |p  in  S5. 

DO  loop  92  stores  the  synthesized  pattern 

N 

[T)f  -  I  o  [Tj<i 
1-1 

in  TF  and  its  magnitude  in  TFM.  DO  loop  92  also  accumulates  ||  [T]f||^  and  the 

pattein  synthesis  error  e  of  (2-23)  in  TFS  and  SI  respectively.  Statement 

* 

56  chtnges  the  angle  of  G  to  that  of  (Tf)  .  At  the  end  of  DO  loop  19,  the 
source  norm  squared 
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the  relative  pattern  synthesis  error  — ,  and  Q  defined  by  (2-15)  as 

WlJ 

M  are  stored  In  FF,  E,  and  Q  respectively. 

l|[Tll|P 

DO  loop  9A  stores  f  in  F  and  accumulates  the  source  norm  squared 

iifip  - !  if„p 

n“l 

In  S3.  DO  loop  94  dlsvovers  that  Che  element  of  ?  Is  largest  In 
magnitude.  As  expected,  the  second  and  fourth  quantities  S3  and  S4 
printed  by  statement  57  are  very  close  to  the  final  values  of  FF  and  Q. 
Statement  59  punches  f  on  cards  for  possible  Input  Into  the  program  of 
section  IX  of  this  report.  DO  loop  Su  normalizes  the  element  of  1 
largest  In  magnitude  to  1  and  then  stores  the  magnitudes  and  phases  in 
degrees  of  the  elements  of  f  in  FM  and  FP  respectively. 

According  to  the  printed  outputs  of  the  magnitude  pattern  synthesis 
with  constrained  source  norm  program  of  this  section  and  the  magnitude 
pattern  synthesis  program  of  section  III  of  this  report,  it  has  been  pos¬ 
sible  to  reduce  the  source  norm  squared  from  27.540  to  about  8.0  with 
only  a  slight  Increase  In  the  pattern  synthesis  error. 
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LISTING  OF  MAGNITUDt  FATTFKN  SYNTFIFSIS  WITH  CONSTRAINFn  SOURCF  NORM  PROGRAH 


//  lOO^A  .FF  ,  ISS,  1  ,4  1  ,  'MAIIT?  .JOE  '  ,REr,inN  =  ?OOK 

//  FXEC  WfTFIV 

//0O.FT02F001  I)D  S  Y  SOOT  =  R,  r)CR=  ( RF  CFM  =  F  ,  RLKS  I  7  E  =  Rn  I 
//Gfl.SYSIN  01)  ’F 

SJOH  MAUT7.  ,T  IMF  =  1  ,PAr,ES=40 

COFiFL  tX  U.Ol  ,CON.I('.,P(  100)  ,T  (  101  ,TP  I  3A0  I  ,0(34) 

COMf LFX  r ( 10) ,ALP( 10) ,TF 134) ,F( 10) 

ni  M(  NSION  X  (  1  0  )  ,Y(  10)  ,PH1  (34)  ,4M(  10)  ,AM2M0)  ,H(  34  )  ,FF(  41  )  ,F  (  41  ) 
f)  I  WE  NS  I  ON  0(  41  1  ,r.f.  (  1  0)  ,  AS(  31  )  ,FS1  30)  ,TFM134)  ,FM(  10)  ,FP(  10) 

PI=: .141443 
CP=f I /IHO. 

0= (C  .  ,  1  .  1 

RFAIll.RO)  N,M,Nr,,N4,RK,FPS 
30  FlIKF  AT  (  4  I  3.2E14 .7  ) 

WRrF!,3,3?)  N,M,N().N4,Rk,EPS 

32  F(IR,FAT('0  N  M  no  nr  '  ,4X  , 'BK*  ,  1  1  X,  'EPS* /1X,4I  3,  2F14.7) 

RFAf'l  1  , 10  1  (  X  (  I  )  ,  I  =1  ,N  ) 

1  0  FOR) A  T  (  1 OFH  .4 ) 

WRP  F(  3,13)  (X(  1  I  .1=1  ,N) 

I  3  FflRNAT  (  '  OX  •  /  (  IX  ,10FH.4  )  ) 

R  FA(l(  1  , 10  H  Y  1  I  1  .  I  =  1  ,N  1 
WR I  ■  F ( 3 , 37 ) (Y ( I  ),I  =  1,N) 

37  FORMAT (  'OY '/( IX, 10FB.4 ) ) 

RFAfM  1  , 101  (PHI  (  I  )  ,I =1 ,M  ) 

WRT  F ( 3.  I  1  I  (PHI  (  I  I  ,I =1  ,M) 

II  FORMAT!  'OPHI '/(IX.IOFH. 4)) 

'■•CAIK  1  ,23  )  (  AM  ( .1)  ,.1  =  1  ,N) 

23  FOPFiAT  (  4F14,7  ) 

WRI'  F  (  3 ,24)  (  AM(,I)  ,  J=1  ,N) 

24  FORMAT (  'OAM'/ ( 1 X,4F 14,7) ) 

NN=M*N 

;m=m 

JAM:  (N*l  )/2 

REAM! 1 ,23) (P( J ) , J  =  I  ,NN) 

WRT  F ( 3,21  )  (P( J) , J  =  1  .NN) 

21  FORMAT!  '  OP  ' / ( 1 X , 4F 1 4 . 7 ) ) 

00  ;r  J=1,N 
AM2 ' J ) = AMI J )4AM ( J I 
24  CON  I  NOE 
on  '4  K  =  1 ,  M 
S  1=)'F  I  (K  )*CP 
CS=(  OS!  SI  ) 

SN=MN(S1  ) 
no  4  I  =  1  ,N 

S  1=|1K<=  (  X  I  I  )*r.S  +  Y  (  I  )«SN  ) 

T(i  =cns(si  i+ofsiNi SI ) 

15  CON  INOE 
Jl  =  ll 
J2=l, 

on  4  J=1,N 
TPl 12 ) =0. 

1)0  7  I  =  1  ,  N 

J1=J1+1 

TP(  j2  )  =  TP(  J2)  +  TI  I  )<=P(  J1  ) 

17  CON  INOE 
J2  =  .I2  +  M 
14  CON  INOE 
14  CON  INOF 
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0(1  22  .)G=1,NG 
PeAD(l,?3l  CN 
WRITE(3,?H)  CN 
28  FORMAT ( 'OCN= ' ,F1  A.T) 
RFADd.inXGd  I, 1  =  1, Ml 
MR  I  TF  I  3,12  I  ir.l  I  I  ,  I  =  1,M) 

12  FORMAT  (  'Of.' /  (  IX,  lOFfl.4  n 
eP=FRS*CN 

r,o=o. 

on  25  .l  =  l,M 
oi=coN.ir,(r,(  j  1 1 
.S1=G(JI»U1 
GG=GG+S1 
H( JI=S0RT(S1) 

G(ji=in 
25  CONTINUF 
FF (11=0. 

F( 1  1  =  1  . 

(T(  11=0. 

JR=1 

AL  =  AM(  .1AM  I 
AS ( 1 )=AL 

201  DO  1<J  KR=1  ,N<J 

.11  =0 

DO  18  .1  =  1, N 

01  =0, 

on  2  0  1  =  1 ,  M 

.il  =  Jl  +  l 

U1=U1  +  TH(.I1  l#G(  1  I 
20  CONT INOt 

C(J)=C0NJG(U1  ) 

cc  (.n=oi*c(  ji 

18  CONTINUF 
S5=CN 

on  31  K=1,N 

SiS  =  S6-CC(K)/AM2(K) 

31  CONTINUF 

IF(S6|  66,85,65 

65  J<»=1 
AL=0. 

AS( J9)=AL 
FS ( JR )=S6 
GO  TO  87 

66  AL=AS(JR) 

AS( 1  )=AL 
JR  =  1 

DO  68  J  =  1  ,30 

S2=CN 

S5  =  0. 

00  67  K=1,N 

S3=AM(K )+AL 

S4=S3*S3 

S2  =  S2-r.C(K  )/S6 

S5=S5+CC(K)/(S4*S3) 

67  CONTINUF 
FSIJR)=S2 

IF(AHS(FS(J9n.LE.FR)  GO  TO  87 
60  AL  =  AS(.IR)-.5*FS(.|RI/S5 
JR=JR+1 

IFIAL.LF.O.I  AL=0. 
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&S(J‘))=AL 
hH  CONTINlJb 

WR I T6 ( I 

21R  FDRMATf  'ONEWTCIN"  MFTHdD  DIO  NOT  CONVERGE*) 

87  S5=0. 

on  <)1  K=1,N 
ALP(K)=C(K)/)AM(K)+AL) 

51  =  ALP  (  K  )*r.CIN.)G(  ALP  (K  )  ) 

SSrSS+Sl 

Rl  CONTINUE 

WRI TE ( 3,90)  ( AS(K ) ,K  =  1 , JR ) 

90  EORMATI*  CONVERGENCE  DE  LAGRANGE  MULT  I  PL  I ER ' / ( 1 X , lOE 1 1 . 9 ) ) 
WRITE(3,100)(ES(K),K=1,J9) 

100  FORMAT!*  FUNCT10N*/(lX,10En.4)  ) 

Si=n. 

TFS=0. 
no  9  2  J  =  1  ,M 
J2  =  J 
U1=0. 

00  93  K=1,N 
U1=U1+TP( J2)»ALP(K ) 

J2= J2+M 
93  CONTINUE 
TF( J)=U1 
lil=CONJG(Ul) 

S2=TF(J)*U1 
TFS=TFS+S2 
TFM( j  )  =SORT (S2 ) 

S2=TFM(J)-H(J) 

Sl=Sl*S2*S2 

36  G) J)=(H(J)/TFM(J) )*U1 

92  CONTINUE 
J2=K9+1 
FF( J2)=S3 

F ( J2  )=S1/GG 
0( J2 ) =7M*S3/TFS 
I F (E ( J2  )  .GE.E(K9 ) )  GO  TO  44 
19  CONTINUE 

44  WRITE(3,5I ) (FF( I ),I=1.J2) 

31  FORMAT (  *ONnRM  SOUAREOOF  F * / (1 X , 3 F 1 4 .7 )  ) 

WRI TEI 3,52) IE( I  I ,  I  =1 ,J?) 

32  FORMAT! *0RFLATIVF  ERROR  * / ( 1 X , 3E 14 . 7  )  ) 

WRITEI3,53)!OII ).I=1.J2) 

33  FORMAT!  *00'/! 1X,3F14.7) ) 

52  =  0. 

53  =  0. 

DO  94  J=1,N 
F! J)=0. 

JI=J 

00  93  1=1, N 
F ! J ) =F  !  J)+P! J1 )*ALP( I ) 

J1 =J1 +N 

93  CONTINUE 
SI=E!J)>»C0NJG!F!J)) 

S3=S3+S1 

IF ! SI  .LT.S2)  GO  TO  94 
S2  =  S1 
J3  =  J 

9^  CONTINUE 
Ul=l./F! J3) 
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S4  =  ?m»c,3/THS 

•>7  W«ITK^,54)  Gr..S^,(-(,l?l,S4 

bi.  FDHMAT  ( '0G(;=  '  ,F  I'* .  7,  •  FF  =  '  ,  F .  7  .  •  F=*,F14.7,'  0=',F1A.7) 

W« 1 TF ( ) (F ( I ) , I =1 ,NI 
4fl  format (  'OF •/ ( 1 X , SF  1A.7  M 
5R  WR I TF (?,?3 ) (F ( n , I =1 ,N) 

on  5b  .1  =  1. N 

F(.l)=Ul*F(.n 
F  M(  J  )=CARS  (F  (.11) 

S1=RFAL(F(.I)  ) 

S2  =  A I  mag ( F ( J  )  ) 

FP(  J(  =  (ATAN?(S?.S1  )  )/f,P 
55  CONT  IMIIF 

WR  I  TF  (3,4?  )  (FM(  J  I  ,.l  =  !  ,N  I 

42  FORMA T (  'OFM> / ( I  X , 1  OF  7.3)  I 
WR  I  TF  (  3,43  )  (FR  (  .1  )  ,.M1  ,N  I 

43  forma  I  (  'off ' / (  1 X  .  1  of  7. 1  I  ) 

WR  I  TF ( 3,44  I  ( TF  (  K  )  ,K  rl  ,M  ) 

4q  FOR'-'A  r  ( 'OTF  '  /  (  1  X  ,  IOFH.4  )  ) 

WR I TF ( 3,50 )  ( TFM (  !  )  ,  I ^ 1 ,M| 

50  FIlRMAl  (  'OTFM  •  /  (  1  X  ,  lOFfl.4  I  I 
22  CONTIMOF 
STOP 
FNO 

tOATA 

10  34  1  40  0,?40H8R?F*00  0. 1000000F-0? 

0,0000  5.340H  10.604H  1‘5.50?7  14.244R  lQ.?f>4«  15.5027  10.404«  5.3iS0R  0.0000 

-10,0000  -4.43?0  -8.4704  -4.?^44  -2.555P  2.5554  4.?844  8,4704  a.43?0  10.0000 

5.0000  15.0000  25.0000  35.0000  45.0000  55.0000  45.0000  75.0000  85.0000  45.0000 

105.00001 15.0000125.00001 35.0000145.0000155.0000145.00001 75.0000185.0000145,0000 
205. 00002  15. 0000225. 0000235. 0000245. 0000255. 0000245. 0000275. 0000285. 000020  5. 0000 
305. 0000 3 15. 0000 3 25. 00003 35. 000034 5. 0000355.0000 
0.4447145F  0?  0.8078324F  0?  0.4145313F  02  0.544P447F  02  0.?B754Q?F  02 

0.2551213F  0?  0.4434283F  01  0.3275412F  01  0.5507557F  00  0.14284236  00 

0.24035226  00  0.2  3044216-04  0 .4844H44F-0 1  0.4 24 1 Q7?6 -08-0 . 344742  IF  00 

-0,47247846-07-0.488551  IF  00-0. ?473 702F-07-0 .  1 98 2043F  00  0. 4377544F-04 
0,19820766  00  0.25471256-04  0.4885510F  00-0. 304 1 387F-04  0.3447433F  00 
-0.1850408F-04-0. 48448846-01  0.27135476-04-0.29035406  00  0.54295756-04 
-0.84149456-01  0. 21413146-04  0.27209406  00-0 . 1  i30?0F-04  0.33087876  00 

-0.44514516-07-0.11724546  00  0. 38587 1 56-04-0. 54 ?4 34? 6  00  0.24475486-04 
-0.54243786  00-0.53487206-07-0.11724736  00-0.52443186-04  0.33087806  00 
0.30878536-04  0.?7?04A1F  00  0.73044446-04-0.89149536-01  0. 2 1 002476-04 
0.52485156  00  0.42702486-04  0.43080886  00  0. 1 4 7 3407 F -05  0.?4?05?36-01 
0.45195296-04-0.77841446-01-0. 55244186-04  0.14351426  00-0.88454346-04 
0.14351506  00  0.1771 575F-04-0. 77842046-01  0.44541546-04  0.24213556-01 
-0.44278556-04  0. 43082036  00-0.14294486-05  0.524B544F  00-0. 7 5? 1 597F-04 
-0.1047320F  00-0.1435871F-05-0.5457034F  00-0.11543896-05-0.39484486  00 
0.14124106-04  0.88930496-01  0 . 44 75079F -04  0.1455435F  00-0. 72441 48F-04 
-0.14554456  00-0. 81522546-04-0. 88929186-01  0. 5?9 ?844F -04  0.34484406  00 
0.48031496-04  0.54549576  00-0.58201976-04  0.  1047221F  00-0.  1 504 1 32F-05 
-0.377347  7  6-01-0.2  7940  796-05  0.18755046-01-0.4  04  ^4046-04-0.41474846  00 
0,34249276-05-0.55731316  00  0.23544546-05-0.11802406  00-0.14773336-05 
-0.11801706  00  0.32594736-05-0.55730485  00  0.55741246-05-0.41477226  00 
0.22424836-05  0.18/58986-01  0. ?4 1390 76 -04-0 . 3 7 7 2 3?8 6-0 1  0.32489496-05 
-0.53597126  00-0.34488546-05-0.11951076  00-0.85048476-04  0.10409196  00 
-0.19182216-05-0.28100816  00-0.51093436-05-0. 32897386  00-0.27118496-05 
0.32897446  on  0.14841946-05  0.28102016  00-0.13471346-05-0.10408116  00 
-0. 32871226-05  0.1  195in?F  00  0.451727  86-04  0.53594926  00  0. 3 1 48825F-05 
-0.41397486  00-0.42519496-04  0.31210946  00  0.35435446-04  0.20240036  00 
0.14254786-04-0.1954)056  00-0.47427016-04  0.38494776  00  0.49871336-04 
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nn-n is‘>407‘»f  no  o.‘)?9444iF-n7  o.?o?40??F  no 

O.;'0‘5iS4n3E-04  0.31?104?F  00  0. 1403'J(SHF-0/S-0.413Q754F  00-0.  ?(SR43?SF-04 
0.32104H3F  no  0.4147647F-04-0.344414?F  00-0 . 4 54 1 9 b 3F -04  0.2344244E  00 
0.22‘37y97F-04  0.1474n03F  00-0 . 494995  1 F-04-0. 4474  349  f  00- 0. 4 1  4fl04  4F-0  7 
().44744n4F  nO-0.4944?04F-n4-0.l475034F  00-0 . 32H R4 54F -04 -0 . 2344 209E  00 
().2H7FH05F-04  n.3444?09F  00-0. 5494532F-04-0. 321 1 029E  00  0 . 59 1 7425E-04 
0.1942254E  on  0.1733023F-05-0.3774?l?F  00-0 . 3023024F -05  0.41S3334F  00 
().3710524F-05-0.3425054F  00-0.  31  2  1  309F-05  0.in  7757F  00  0.  1  425577F-05 
0.11  17430F  00-0.4?75414F-04-0.34249?5F  00-0  .  «  2  7H  H  49F -07  0.41»'3224E  00 
0.iy3R445F-0h-n.3774149F  00-0. 35R4534F-04  0.1942231F  00  0. 1 24 1 534F-07 
-0.119H0I5F  00-0.  455H49K-04  0.2539174F  00  0 . 4 1 11 4 M 7F -04 -0 . 3 79 1 949E  00 
-0. 3514973F-04  0.39n7340F  00-0. 3590495F-04-0. 35311 59F  00  0. R05373 1 E-04 
0.3531210F  00-0.-  3fl4449e-05-0.39n74«4F  00  0 . 294 R 4 2 2 F -05  0.37920«4F  00 
-0.342451 1E-05-0.2539?79F  00  0. 2444993F-05  0.119R04RF  00-0. 1 144047F-05 
0  .MOOOOOOF  +  01 


4.0191 

n.onon 

3.443  / 

0.0000 

2.3442 

0.0000 

1 .7434 

0.0000 

1.414? 

o.oono 

1  .220H 

o.nooo 

1 .1034 

0.0000 

1 .0353 

o.nooo 

1 .0034 

0.0000 

o.oono 

0.0000 

o.nooo 

0.0000 

0.0000 

0.0000 

0.0000 

o.nooo 

o.oono 

0.0000 

0.0000 

o.nooo 

0 .0000 

o.nooo 

o.nooo 

0.0000 

0.0000 

o.nooo 

0.0000 

0.0000 

0.0000 

o.oono 

o.nooo 

o.nooo 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.nonn 

o.oono 

o.nooo 

o.oono 

o.onon 

0.0000 

o.nooo 

0.0000 

0.0000 

o.oono 

o.onon 

0.0000 

0.0  00 

o.oono 

0.0000 

0.0000 

0.0000 

0.0000 

o.oono 

0.0000 

o.oono 

('  .0000 
iS  TUP 
/• 

// 

0.0000 

PRlNTET  OUTPUT 

N  M  NO 

N9 

BK 

EPS 

10  34  1 

40  r..?90HHH?p+00  O.lOOOOOOF-02 

X 

0.0000 

5.3404 

10.4044 

15.5027 

19.2444 

19.2444 

15.5027 

10.4044 

5.3404 

o.nooo 

Y 

-10.0000 

-9.4'-20 

-4  .4704 

-4.2494 

-2.5559 

2.5559 

4.2404 

4.4704 

9.4320 

10.0000 

PH  I 

5.0000  15.0r00  25.0000  34.0000  45.0000  55.0000  45.0000  74.0000  45.0000  95.0000 
105.00001 15. 0(00 12 5. 0000 135. 00001 45. 0000 155. 0000 14 5, 0000 175. 0000 145.0000 195. 0000 
2U5.00002 1 5.0(00225.0000235.0000245.0000255.0000245.0000275.0000244.0000295.0000 
305.n000315.n(  0032  5.0000  33  5.0000345.0000355.0000 

AM 

n.9447145F+0;  0.H07H329F+02  0. 4 1953 1 3F+0?  0 . 544449 7 E +02  0. 247599?F+02 
0.255121 3F  +  0:  0.94392H3F  +  01  0 .32754 1 2F  +  ni  0.5507557E+00  0. 142«423E+00 


0.2903522F  +  0(  0 . 23044? 1 F-04  0. 4R44444F-0 1  0 . 4?4 1 9 72E-04-0. 344 742 1 F +  00 
PLUS  39  MOHE  I  INES  OF  PRINTED  OUTPUT 

r.N-  0.40000001  +01 
(, 


5.0191 

0.0(00 

3.4437 

0.0000 

2.344? 

0.0000 

1.7434 

0.0000 

1 .414? 

o.nnof 

1 .220H 

o.onon 

1 .1034 

0.0000 

1 .0353 

0.0000 

1 .0034 

0.0000 

0.0000 

n.noo' 

0.0000 

o.nooo 

0.0000 

o.oono 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

o.onni 

o.nooo 

o.nooo 

0.0000 

o.oono 

0.0000 

0.0000 

0.0000 

0.0000 

o.oono 

O.OOOi 

0.0000 

o.oono 

0.0000 

n.onoo 

0.0000 

0.0000 

o.onon 

0.0000 

0.0000 

0.000' 
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o.onon  o.ooon  o.ooon  o.oooo  n.onoo  o.oooo  o.oono  o.oooo  o.nooo  o.noon 

O.OOOn  0.0000  0.0000  o.oooo  o.oooo  o.oooo  0.0000  0.0000  0.0000  0.0000 

0.0000  o.oooo 

CDNVt-ROt-NCl-  1)1-  LAOkANlit-  M1JLT)'LI(:K 

o.?H7^(-+o2  n.onnoh  +  no  o.<.OHOh-oi  o.ft?3?e-oi  o.isisist-ri 

HtINi:  1  ION 

0.77  ^lE+Ol-O.bHAbh+Ol-O.  14  7/, F +01-0.  l‘<77F+00-0.47596-02 
CONVI-ROFNCt  OF  LA(.KAN(;f  MULTPLIFK 
O.hhlRF-Ol  n.l014F+00  n.lU4F+00  0.11?1F+00 
FONC I  ION 

-0.22FlF  +  01-0.4  1HbF  +  00-n.?n?F-01  -0.IS294F-04 
CONVFHOFNCF  OF  LAOHANGF  MOI.TPLIFK 
0.1121F+00  O.lPbbF+OO  O.1P4SF+00 
Forjc  T  ION 

-0.44  31 F  +  00-0.3b41  F-ni  -n.  1  7  0  7F-03 
CONVFHGFNCF  OF  LAGRAN)))-  MOLTFLIFK 
0.124bF,00  0.1734F+n0 
FONC  T  iriN 

O.H7‘)3F-01-0.1  lObt-02 
CONVERGENCE  OF  LAGRANGE  MOl.TPLiFR 
0.1234E+00  O.llbflE+00  n.ll62F+00 
FUNCTION 

O.?7n‘)E+0n-0. 1  127E-01-0.  lHlPF-04 
convergence  OE  LAGRANGE  MOLTPLiFR 
0.11h2F+00  n.l?73E+00  0.12HOF+00 
EONCTirjN 

-0.42H41 +nO-0.?3|hE-01  -n.7h?SF-04 
CONVEKGENM  HE  LA(,RAN(.E  M|)|  1PlI6« 

o.i2RnE,iin  fi.is?nF  +  oo  o.i4S7F+oo 
FONC  T lOM 

-0.4443)  ♦0li-(  .  I  OOH  F  +  nO-O.  14  1  OE-02 
CONVERGENCE  OF  LAGRANGE  IAIJLTPLIER 
0.14‘)2E+(ili  O.147?F  +  n0  0.1A)7HF+00 
F)iNCTION 

-().347bF/0O-0.7()0AiF-01  -n.SHl  7F-04 
CONVERGENCE  OE  LAGRANGE  MOLTPLiFR 
0.147RF+n0  0.1704F+00 
FONC T ION 

-0.434‘)F-01-0.1?17E-n? 

CONVERGENCE  OE  lAGRAO.F  MOLTPLIFR 
0.1704E+00  ().1713E  +  n0 
FONCIION 

-0.1031E-01 -0. i407e-n4 
CONVERGENCE  OE  LAGRANGE  MOLTPLIFR 
0.1713E+00  0.17n4F+n0 
FUNCTION 

0.1030E-01-0.1421F-04 
CONVERGENCE  OE  LAGRANGE  M))LTPLIFR 
0.1704E+00  0.17nbE+n0 
FONCT  It)N 

0.11b4E-01-0.?003F-04 
CONVERGENCE  OF  LAGRANI.E  M))LTPLIFR 
0.170bE+00  0.1702E+00 
FONCT ION 

0.RRt)2E-02-0.1  144E-04 
CONVERGtNCE  OE  LAGRANGE  MOLTPLIFR 
0.1702E+0n 
FONCT I  ON 
0.44  /4E-02 
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MOKM  SOIIAHU)  (Ih  t- 

0 .nnnoonot-  +  nn  t) .flonATb?^  +  oi 
o.Hnnoni3t+f)i  o.HnnoohHt+(n 
().HnnoooRi-+oi  n.noonno'jH +ni 


o.HooooRAF+ni 
O.HO014O?(-*Ol 
O.H00001 lh+01 


o.Rnoni63i-+ni 
o.Ronnnsoh+ni 
o.«nnor)07t'  +  ni 


0.  Honio'iRt^  +  oi 
n.«noi?OHF+oi 
0.7qq40?lR+01 


Wpl.ATIVt- 

o.iononooH*n) 

().?3301  h2F+0r 
0.141  +nf 


o.?40H044F+nn 

o.?\io?onF+on 

0.1411 200F+00 


n.2hAR3R4F'*on 
0.  14S4«''6F  +  00 
0.14111?7F+0n 


0.?R4?nH4,F  +00 
0. 14?ni 7hF+00 
0.1411 1041+00 


O.^ARFIOOF+OO 
0. 141?«03F+00 
0. 191 1?H^F  +00 


O.OOOOOOOF+0( 
0.  n  90:)70F  +  01 
O.hH^hFRSF+Ol 


0.77H7 /OHh+ni 
0.fe44  74H  7F  +  01 
0.feH?‘)104F+01 


0. 74455??F+01 
0.687Hft3?F+01 
0.6H?3477F+01 


0.737>34‘3SF  +  ni 
0.M<3‘)937F+01 
0.7>«?31  ^7F  +  01 


0.7?R7‘iiSRF+01 

0.lSR?9R?4F+0l 

0.ftH17fl93F+01 


0(;=  0.‘5‘>S9404t  +0?  FF=  0.79434446  +  01  F=  0.19n?H4F+n0  Or  0 . 4H 1  7R?9F +  0 1 


F 

-0.4l7731Mt  00-0.4?49Rb41--01  0.H441H43F  00  0.4310747.6  00-0 . 94 1  ?030F  00 

-(I.110O147F  O' -0.96H7  7476-01  0.1414407F  01  0.4497,797,6  00-0. 4?93?H7,F  00 

-0.440R440F  00  0.477,77,076  00  0.310H7,316  00-0 . 44  H  6  H9  76  00 -0  .  RR9 1 094  6 -0  1 

0,H1724?16  00-0.330941  06  00-0.777,40346  00  0.47,4701  46  00  0.17R90776  00 


6M 

0.240 

0.64,  1 

1.000  0.4H1  0.427  0.424  0.377 

0.468 

0.44? 

0.333 

6  6 

-93.7, 

147.  > 

-0.0  -1 34.4  86 

.6  -43. 

4  74.7 

-133.7 

16.1 

146.1 

T6 

? .81 84 

0.1  76H 

3.2367  -0.3112 

7.4743 

-1  .0900 

0.096? 

-1  .4747 

-0.6606 

-0.904 

-1 . 6440 

0.36H4 

-0.784?  1.3041 

0.7037 

1 .7408 

0.6494 

0. 5049 

0.4381 

-0.061 

O.lOl 1 

-0.1 404 

-0.03HH  -0.0344 

-0.0711 

0.0096 

-0.1 747 

-0.0149 

-0.1940 

-0.062 

-0.2238 

-0.114? 

-0.147?  -0.143R 

-0.1 138 

-0. 1 183 

-0.0407 

-0.0468 

-0.0764 

-0.074 

-0 .07,8  1 

-0.0637 

-0.1 186  -n.l4?4 

-0.1447 

-0.200? 

-0.1779 

-0.7330 

-0.1 768 

-0.748 

-0.01 R3 

-0.1 366 

0.0767  -0.0381 

-0.046? 

0.0483 

-0. 1777 

-0.0039 

-0.106? 

-0.076 

-0.17,2  3 

O.lORl 

-0.41(1?  0.7446 

-0.6066 

0.7081 

-0.3449 

0.0903 

0.4943 

0.108 

1.7141 

0.7472 

TFX 

2.8247 

3.2416 

7.800?  1.7747 

1.170? 

1.  919 

1.4764 

1.7474 

n.8??4 

0.44? 

0.1812 

0.0446 

0.0718  0.1764 

0.2040 

0.7417 

0.7401 

0.1641 

0.0699 

0.03'= 

0.0932 

0.1844 

0.7437  0.7401 

0.2741 

0.1874 

0.0466 

0.0741 

0.1278 

0.109 

0.1440 

0.4781 

0.6413  0.366? 

0.4071 

1.7338 
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VII.  PATTERN  SYNTUKSIS  WITH  CONSTRAINKD  QUALITY  FACTOR 


The  program  of  this  section  Is  obtained  by  replacing  the  block 
statements  from  (and  including)  statement  201  to  statement  91  in  the 
pattern  synthesis  with  constrained  source  norm  program  of  section  V  of 
this  report  by  the  block  of  statements  listed  in  the  present  section. 

The  punched  card  input  data  is  the  same  as  that  of  the  pattern  synthesis 
with  constrained  source  norm  program  except  that  CN  is  interpreted  not  as 
the  upper  bound  on  the  source  norm  squared  but  as  the  upper  bound  on  the 
quality  factor  Q. 

Summarizing  the  results  of  section  VI,  Report  No.  2,  the  optimum 
source  vector  f  is  given  by 


Q  N  C, 

►  «  1  ► 

fa  (1  4*  £  — )  )  ■  t) 

“  iil  'i  1 


(2) 


where  Q  is  the  upper  bound  on  0,  C,  is  given  by  (2-42)  and  A,  and 

O  1  11 


satisfy  (2-35).  If 


0 

o 


then  6*0.  Otherwise,  8  is  determined  from 


(3) 


If 


N 


i.i 

^o  -  ”  N'lcj^A-;  =  =  ° 

i^i 


(4) 


(5) 
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where  Is  the  largest  eigenvalue  of  (2-35)  then  (4)  has  no  solution 
because  Q  Is  always  greater  than  Q^.  If 


•^1  1  <^0  ^  «  N 


N 

I  |cj- 

1-1  ^ 


I 

i-1 


(6) 


then 


-a>  <  B  £  - 


If 


'^2 


0 

o 


(7) 

(8) 


0  •  6  ■  » 


(9) 


The  above  scheme  determines  B  uniquely  because  the  left  hand  side  of 
(4)  is  a  monotone  decreasing  function  of  8  in  the  intervals 
and  (0  “) , 

a  verbal  flow  chart  of  the  replacement  block  of  statements  is  as 
follow;!.  Just  after  DO  loop  202,  and  of  (6)  and  (3)  are  stored 

in  Q2  and  Q3  respectively.  Starting  with  either  8  =  -  A^^  - 
3  -  ^(■j+l)/2’  loop  213  carries  out  Newton's  method  for  finding  the 
root  o'  (4).  According  to  Newton's  method. 


F(S  ) 

®I+1  °  ®1  ■  F'{6j) 


(10) 


differentiating  (4),  we  have 


F'(6) 


2M 


N  N  IcJ^A^  N  N  |cJ2 

1=1  (^1+3)^  i-iV^ 


! _ , 

i-i 


(11) 
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i 


An  alternate  form  of  expreaaion  (11)  Is 


F'(6) 


2M 


1^1  yip  "  1-1 


N  |c^ 

1-1 


(12) 


For  large  |  g|,  (11)  Is  the  difference  of  two  terms  each  of  which  Is  proportional 
to  -T,  but  (12)  Is  the  difference  of  two  terms  each  of  which  Is  proportional 
to  Thus  (12)  Is  less  susceptible  to  round  off  error  when  |s|  Is  large. 

P 

However,  both  (11)  and  (12)  are  Indeterminate  when  3  -  -  Equation  (12) 

can  be  written  as 


F’(6) 


2H  (,  1^1’  ,  ? 

/  ?  1  1\2 

lii 


MA^+B)^  (A^+6)^^  ‘  MA^+S)^ 


(13) 


After  cancellation  of  the  terms,  expression  (13)  becomes 

tA 


F'(R) 


2M 


^  (a'  +6)  ^  ^  7ir+«73  A  L  (X  +r^jA 


1-2 


1-2 


(A^+B)^  ^^2  ( V®’ '  ^  ill  1  i.2 


N  ICJ^A 


51 


DO  loop  214  accumulates  in  SI,  S2,  S?,  S4,  and  S5  the  five  different 
sums  appearing  In  (4)  and  (14).  Next  F(B)  of  (4)  is  stored  in  FS  and  F'(B) 
of  (14)  In  DO.  B  of  (10)  is  stored  in  both  Al.  and  AS.  If  the  new  B  is  out¬ 
side  the  correct  interval,  the  logic  between  statements  216  and  218  replaces 
B  by  the  average  of  the  previous  B  and  the  end  point  of  the  interval.  DO 
loop  91  stores  the  coefficients  of  ip^  appearing  in  (2)  in  Al.P. 

The  sample  punched  card  input  data  is  numerically  Che  same  as  that 

of  the  pattern  synthesis  with  constrained  source  norm  program  except  that 

CN  =  C . 4000000  E+01  is  replaced  by  CN  =  0.4049586  E+01.  The  number 

0.404S586  E+01  obtained  from  the  printed  output  of  the  pattern  synthesis 

* 

with  constrained  source  norm  program  is  the  Q  of  the  source  vector  f  which 
minimizes  the  pattern  synthesis  error  subject  to  the  constraint  ||f||^^  4. 
Comparing  printed  outputs,  we  see  that  the  pattern  synthesis  error  of  the 
pattern  synthesis  with  constrained  Q  program  is  only  slightly  smaller  than 
that  cf  the  pattern  synthesis  with  constrained  source  norm  program. 
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Tf)  flHT/klN  THE  PATTERN  SYNTHESIS  WITH  CONSTRAINEf)  0  PROGRAM,  REPLACE 
THE  HLnCK  C)E  STATEMENTS  FROM  STATEMENT  201  TTI  STATEMENT  91  IN  THE 
PATTERN  SYNTHESIS  WITH  CONSTRAINED  SdllRCE  NORM  PROORAM  RV  THE 
EOLLnwlNG  STATEMENTS. 


;M2=2.*2M 
AMM=-AM( 1 ) 

01 =2M/AM( 1 ) 

S1=0. 

S2*0. 

53  =  0. 

54  =  0. 

on  20?  K=1.N 
S1=S1+CCIK l/AM?(K ) 

S?=S2+CC(K)/AM(K) 

S3=S3+CC(K ) 

S4=S4+CC(K I^AMIK ) 

202  CONTINUE 
02=ZM*S3/S4 
03=ZM*S1/S? 

WRITE(3,?03)  01.02,03 

203  FORMAT! '001*' ,E14. 7, •  02=*,F14.7,'  03=',E14.7) 
IFICN.OE.Ol )  GO  TO  204 

WRlTE(3.?0b) 

205  FORMAT! 'OSPECIFIED  0  IS  TOD  SMALL') 

GO  TO  2? 

204  J9=l 

IF!CN.LE.03)  GO  TO  207 
AL=0. 

AS!  J9)=AL 
FS!  J9)=CN-03 
GO  TO  POR 

207  IF!CN.GT.O?)  GO  TO  211 
K0  =  1 

AL=-AM!1)-AM!JAM) 

GO  TO  21? 

211  K0=2 
AL<:AM!  JAM) 

212  AS!J9)=AL 

on  213  K=l,30 
S1*0. 

S2=0. 

53  =  0. 

54  =  0. 

S5*0. 

IFIN.EO.l)  GO  TO  225 

on  214  J=?,N 

S(S  =  AM!J)+AL 

S7=S6*S6 

SH»CC! J)/S7 

S9=SR/S4 

Sl=SI+Sfl 

S2*S2+SR*AM! J I 

S3=S3+S9 

S4*S4+S9*AM! J I 

55  =  S5  +  S9>»AM?!J) 

214  CONTINUE 

225  S6=AMI 1 )+AL 
S7=S4*S6 
SH=CC!1 )/S7 
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!5l=Sl+SH 

s?  =  s;'+s«>i‘AM(  1 1 
S‘J  =  S»/Sh 

HS(  jwi=r.N-;M#si/s? 

IF(AHS(^•S(.N(  ).LF.(-P)  (;n  TO  ?0« 

nO  =  ZM?*  (  +  (  Sb+ AM?(  l)‘i‘S^-?.'i‘AM(n*SA))/(S?*S?l 

al=:AL-fs(.|4)/(ho 
oil  TIK  ?1  h,  ?1  V  I  ,K(0 
?  16  I  F  (  AL  .LI  .  AMM  )  OH  Tl)  ?1  H 
AL  =  .*i«(/MM+AS(.l4)  I 
0(1  1(1  ?  H 

I  F  (  AL  .01  .0.  )  0(1  Til  ?1  H 
AL  =  .t>»A‘,(.l4| 

?|H  .|4  =  jq  +  l 

AS(  Ji)  (=,iL 
2\^  CUNTINOI 

(JK  I  TF  (  ^  .?14  ) 

2\^  Fni<MAT  (  '  OMFWTIINS  MFTH(in  H I  f)  NOT  OINVFROF') 

2n«  s^=i  .  +  Ai.i'r,N//M 
on  W1  K:  1  ,N 

ALP(K)=  S^/(  AM(K  l+AL  I  )*r,  (K  ) 

41  CriNTINIK 

pRIMTFI)  mninil  of  PATTFRM  SYNTHFSIS  with  r.nNITRAINFD  1.1  PRnORAI»l 
N  M  MO  PK  FPS 

in  ih  1  o.?‘inHHH?F+nn  o. looonooF-o? 


X 

o.oono  (j.SOOP  10.60AH  IS. 5027  19.264R  1Q.264R  15.5027  10.604F  5.360H  O.OOOf 

Y 

-in. 0000  -4.6(20  -8.4706  -6.2844  -2.5550  2.5550  6.2R04  8.4706  0.6320  10.0001 

PH  I 

5.0000  15.0000  25.0000  35.0000  45.0000  55.0000  65.0000  75.0000  85.0000  45.0001 
105 .00001 15. 0)00125. 00001 35. 000014  5. 0000155. 00001 65. 0000175. 00001 85. 0000 145. 000 
205 .00002 15.0000225.0000235.0000245.0000255.0000265.0000275.0000285.0000245.000' 
305 .0000315.0  100325.0000335.0000345.0000355.0000 

AM 

0.4467145F+n  '  0.H07R324F+02  0 . 6 1453 1 3F  +  02  0. 546H447F+02  0.28754428  +  02 
0.25512136+0?  0.4634283F+01  0 . 3275412F+01  0.5507557F+00  0. 1628423E+00 

P 

0.24035226+00  0. 2304621 F-06  0.68448666-01  0. 424 14726-08-0.  364  742 1 F  +  00 
PLUS  34  MOHF  .INFS  OF  PRINTFO  flllTPIIT 

CM=  0 .4044586 -+01 

0 


5.0141 

0.0  100 

3.8637 

0.0000 

2.3662 

0.0000 

1.7434 

0. 0000 

1 .4142 

0.000 

1  .2208 

0.0)00 

1 .1034 

0.0000 

1 .0353 

0.0000 

1 .0038 

0.0000 

0.0000 

0.000 

0.0000 

0.0  )00 

0.0000 

0.0000 

o.oono 

0.0000 

o.oono 

0.0000 

0.0000 

0.000 

0 .0000 

0.0)00 

0.0000 

0.0000 

o.oono 

0.0000 

0.0000 

0.0000 

0.0000 

0.000 

o.oono 

0.0  )00 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.000 

o.oono 

0.0)00 

0.0000 

o.oono 

0.0000 

0.0000 

0.0000 

0.0000 

o.nnoo 

0.000 

0.0000 

0  .noon 

0.0)00 

0.0)00 

0.0000 

0.0000 

0.0000 

0.0000 

o.oono 

0. 0000 

0.0000 

o.nnn 
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cc 

0.2!jt>2ft61B  +  03  0.4?031bAR+n3  0.?R340ftlF  +  0?  0.  ?70iS  1  2HF+03  0.  3(S1  31 1  5F+03 
0.1311416E+03  0.1  O.lPROlPbF  +  O?  0. 257073QE +00  0. 2735 140F +00 

01=  0.3H02624E  +  00  «?=  0 . h06332(SE  +  00  03=  0 . 1 2t>‘)200E +0? 

CONVEKOFNCE  CIF  LAGRANGF  MULTIPLIER 

0.2H7RE+02  0.143flE+0?  0.7190F+01  0.3‘)V5F+01  0.17Q7E+01  0.RPH7F+00  0.44P4F+n0 
0.107PF+00  0.1878F  +  00  0.23S‘>F  +  00  0.?4<SOF+00 
FUNCTION 

0.3182F+01  0.3068E  +  0I  0.2P13F  +  01  0.2(SR7iF+01  0.2384F  +  0I  O.lRP^iF  +  Ol  O.lllRF+OI 
-0. 2231 E  +  01 -0.6632 E  +  OO-0. in34F +00-0. 3842F-02 

00=  0.533450‘JE  +  02  FF=  0.4??8370E+01  F=  0.323B363F  +  00  0=  0.4032PP0F +0 1 


-0.3523602E  00  0.24637H2E  00  0.8174110F  00  0.1P41064F  00-0. 9640RPRF  00 
-0.33972816  00  0.16696346  00  0.8117120F  00-0.6 356643F -03-0 . 2 1 8fl625F  00 
-0.3119421E  00  0.48239226  00-0.  1 727355F-01-0.  1 R73600F  00-0. 20161  32F  00 
0.6348063E  00-0.37061206  00-0 .444396 1 F  00  0.4400374F  00-0 . 1 763967F  00 


FM 

0.421 

0.822  1 

.000  0. 

811  0. 

214  0.362  0.184 

0.652 

0.  708 

0.4  64 

6P 

-54.4 

17  3.9 

0.0  -121.0  70.4  -76.5  65.3 

1 

■0 

• 

X 

18.5 

13R.7 

TF 

2.4783 

0.0978 

2.9884 

-0.0485 

2.7169  -0.3176 

1.8472 

-0.3935 

0.9372 

-0.0868 

0.5064 

0.3740 

0.371  1 

0.5426 

0.7067  0.2R39 

0.6150 

-0.0921 

0.4150 

-0.2912 

0.3043 

-0.3618 

0.2374 

-0.4543 

O.ORRR  -0.5405 

-0. 1129 

-0.5230 

-0.2497 

-0.4307 

-0.2671 

-0.3355 

-0.1894 

-0.3448 

-0.0762  -0.3690 

-0.0033 

-0.3834 

-0.0332 

-0.3860 

-0.1338 

-0.4174 

-0.2637 

-0.4997 

-0.2536  -0.5H64 

-0.1027 

-0.5931 

0.1114 

-0.4642 

0.2383 

-0.2101 

0.2222 

0.0673 

0.0100  0.2100 

-0.2297 

0.1669 

-0.3980 

0.0737 

-0.5693 

0.0569 

-0.7933 

0.0591 

-0.8963  -0.C296 

-0.5813 

-0.1223 

0. 2582 

-0.1086 

1 .4194 

0.0250 

TFM 

2.4804 

2 .9888 

2.7334 

1.8887 

0.9412  0.6295 

0.7878 

0.7616 

0.6218 

0.5069 

0.4728 

0.5126 

0.5477 

0.5370 

0.4979  0.4446 

0.3935 

0.3768 

0.3834 

0.3874 

0.4448 

0.5650 

0.63H9 

0.6019 

0.4774  0.3330 

0.2322 

0.2102 

0.2839 

0.4048 

0.5721 

0.7977 

0.8963 

0.5941 

0.2801  1.4196 
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Vlll.  MAGNITUDK  HATTKRN  SYNTHKSIS  WITH  CONSTRAINED  QUALITY  FACTOR 


The  program  of  this  section  Is  obtained  by  replacing  the  block  of  state¬ 
ments  from  (and  Including)  statement  201  to  statement  91  in  the  magnitude  pat¬ 
tern  syntf.esis  with  constrained  source  norm  program  of  section  VI  of  this  report 
by  the  bl('Ck  of  statemetns  listed  in  the  present  section.  The  punched  card  Input 
data  is  tie  same  as  that  of  the  magnitude  pattern  synthesis  with  constrained 
source  noi-m  program  except  that  CN  Is  Interpreted  not  as  the  upper  bound  on  the 
source  noim  squared  but  as  the  upper  bound  on  the  quality  factor  Q. 

A  verbal  flow  chart  of  the  replacement  block  of  statements  Is  as  follows. 
Steps  1,  and  3  of  the  iteration  procedure  described  in  section  IV,  Report 
No.  2  are  executed  for  the  K9*'^  time  inside  DO  loop  19.  DO  loop  18  stores  Cj 
of  (2-42)  In  C(J)  and  In  CC(J).  Just  after  DO  loop  202,  and  of 

(6)  and  (1)  are  stored  in  Q2  and  Q3  respectively.  The  logic  between  statements 
203  and  2  .2  adjusts  the  starting  value  of  B  appearing  in  (4)  If  it  is  not  in  the 
correct  interval.  6  is  stored  in  both  AL  and  AS.  For  K9-1,  the  unadjusted 
starting  "alue  of  6  Is  ‘  unadjusted  starting  value  of 

6  is  the  -'alue  of  8  calculated  from  the  scheme  (3)-(9)  for  the  previous  value 
of  K9.  DO  loop  213  carries  out  Newton's  method  for  finding  the  root  of  (4), 

DO  loop  2.4  accumulates  in  SI,  S2,  S3,  .S4,  and  S5  tlie  five  different  sums  appear¬ 
ing  in  (4  I  and  (14).  Next,  F(S)  of  (4)  is  stored  in  FS  and  F'(8)  of  (14)  in  DQ. 
If  the  nev  f  calculated  from  (10)  is  out.slde  the  correct  interval,  the  logic 
between  s  ;atements  216  and  218  repl.aces  fi  by  the  average  of  the  previous  6  and 
the  end  pjint  of  the  Interva'.  DO  loop  91  stores  the  coefficients  of  appear¬ 
ing  in  (2  p  in  AI.P. 

The  simple  punched  card  Input  data  is  numerically  the  s.iroe  as  that  of  the 
magnitude  pattern  synthesis  with  constrained  source  norm  program  except  that 
CN  =  0.8030000  K+01  is  replaced  by  CN  =■  0.6817829  l-Xl! .  The  number  0.6817829 
F.+Ol  obtained  from  the  printed  mil  put  of  the  magi'.itiide  p.ittern  synthesis  with 
constrained  source  norm  progr.am  is  tin'  n  f  the  source  vector  f  which  minimizes 
the  pattern  synthesis  error  subject  to  the  constraint  ^  8.  Comparing 

printed  oitpiits,  the  pattern  synthesis  eri'or  of  the  magnitude  pattern  synthesis 
with  constrained  0  program  is  onlv  slightly  sm.aller  than  that  of  the  magnitude 
pattern  synthesis  with  constrained  source  norm  progiam. 


TO  OBTAIN  IHt  MACNITune  PATTERN  SYNTHESIS  WITH  CONSTRAINED  D  PROGRAM, 
REPLACE  THE  BLOCK  OE  STATEMENTS  EROM  STATEMENT  201  TO  STATEMENT  91 
IN  THE  MAGNITOEU  PATTERN  SYNTHESIS  WITH  CONSTRAINED  SOURCE  NORM 
PROGRAM  BY  THE  FOLLOWING  STATEMENTS. 

AMM=-AM( 1  ) 

7M2=2.*2M 
01  =  ZM/AM(  1  ) 

IE  (CM.fiE  .01  I  GO  TO  20^ 

WRITE(H,20SI 

20S  FORMAT ( 'OSPECIEIEO  O  IS  TOO  SMALL') 

GO  TO  2? 

20A  0(1  19  K9  =  l  ,N9 
J1=0 

no  IH  .1=1, N 

01=0. 

on  20  1=1 ,M 

.)l=Jl  +  l 

01  =ui  +  TP ( .11  )*(;(  1 1 
20  CIlNTINOE 

C(j|=coNjr.  (01 1 
CC(J)=01*C(.I) 

IR  CONTINUE 
SlrO. 

S2=0. 

ST  =  0. 

S4  =  0. 

0(1  202  K  =  1  ,N 

51  =S1+CC(K  )/AM2(K ) 

52  =  S2+Cr.(K  l/AM(K  I 
S3=S3+CC(K I 

SA  =  S4  +  CC ( K  )«4M(K ) 

202  CONTINOE 
02=ZM*S3/S4 
03=ZM*S1 /S2 
WRITE(3,?03)  01.02,03 

203  FORMAT! '001=', E14. 7,'  02=',F14.7,'  03=',E14.7) 

J9=l 

IF (CN.LE  .03  I  GO  TO  207 
AL  =  0. 

AS( J9)=AL 
FS(  J9)=C('J-03 
GO  TO  208 

207  IF(CN.GT.02I  GO  TO  211 
KfJ=l 

IF  (  AL .GT. AMM )  AL =-AM ( 1 )-AM ( J AM  ) 

GO  TO  212 

211  K{J  =  2 

IF(AL.LT.O.I  AL=AM(JAM) 

212  AS(J9)=AL 

00  213  K=l,30 
S1=0. 

S2  =  0. 

S3=0. 

54  =  0. 

55  =  0. 

IF(N.Fn.l  I  GO  TO  22'j 
on  214  .1  =  2, N 
S4  =  AM(  .1  )  +AL 
S7=S4«S6 
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S«=CC(.I  )/S7 
S^sSH/S* 

SI =S1+SH 

S2  =  S2  +  S8*AM( .]  ) 

S4.  =  S44-S‘<>i‘AM(  J  ) 

S5  =  Si)  +  S‘<*AM?(  .)  ) 

21A  CnNTINDf 
22S  S8=AM( 1 ) ♦AL 
S7  =  Sh>i'Sft 
SH=CC  (  n/S7 
SlrSlfSH 
S2  =  S2  +  SH*AM( 1  I 
St>  =  S0/Sft 

FS(  J‘))=CN-ZMWS1/S2 

if{A8S(8S(>‘  ij.le.epi  r.n  in  20« 

|)0  =  ZM?«(  SB'  SS-SA’i'SA+S'J*  (  S5*AM?(  1  )  •S3-2 .  <•  AM  (  1) #SA  I  1  /  (  S2* S2 1 
AL  =  AL-FS  (  J‘* )  /no 
GO  TO  (210,217  )  ,K(J 
210  IF ( AL .L 1 .AMM  )  GO  TO  2l« 

AL*.*)*!  AMM+AS(  J‘<)  1 
GO  TO  21« 

217  IF ( AL .GT .0.  1  GO  TO  21H 
ALs.S'i'ASl  JM) 

21fl  J<J  =  J‘1  +  1 

AS(J<J  I  =  AL 

213  CONTINOe 
W«ITF(3,21V| 

214  FORMAT! 'ONtMTUNS  METHOD  DID  NOT  CONVERGE*) 

20H  S3=l .+AL»CN/ZM 

S5  =  0. 

DO  9  I  K  =  1  ,  N 

ALP(K )=(S3/(AM(K)+AL) I^CIK) 

SA  =  ALP(K  )4rONJG( ALP(K)  I 
SS=SS+SA 
91  CONTINUt 

PRINTED  OUTPUl  OF  MAGNITUDE  PATTERN  SYNTHESIS  WITH  CONSTRAINED  0  PROGRAM 

N  M  NG  N9  BK  EPS 

10  30  1  40  ( .2908Hfl2E+00  0. 1 OOOOOOF-0? 


X 

o.onoo  S,3<08  10.6048  1S.S027  19.2648  IQ. 2048  15.S027  10.6048  5.3608  O.OOOI 


Y 

-10.0000  -4.6-20  -8.4706  -0.2844  -2.5559  2.5559  6. 2894  8.4706  9.6320  lO.OOOt 


PHI 

5.0000  15.0(00  25.0000  35.0000  45.0000  55.0000  65.0000  75.0000  85.0000  95.0001 
105 .0000 115. OC  001 25 .00001 35.0000145.0000155.000016  5.00001  75. 00001 8  5. 00001 95 .000( 
2  05 .00002  15.0(00225. 0000235. 0000245.0000255.0000265. 0000275. 000028  5. 0000295. 0001 
305. 00003 15. or  003  25.0000335.000034  5.0000355.0000 

AM 

O.4467105E  +  O;  0.807H324E  +  02  0 .6 1 953 1 3F +  02  0 . 5468497E +02  0 . 2875992E +02 
0.25512  1 3b  +  0;  0 . 9634  2M 3E +01  0. 32754 1 2E  +  0 1  0. 5507557E+00  0. 1 628423F+00 

P 

O.2403522h+no  n.23O4021F-O6  0 .6844866E-01  0 .4 24  1 972E-08-0 . 364 742 1 E +00 
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PLUS  3*1  MOHF  LINfcS  OF  PK  INTRO  lUITPUT 


CN=  0.6«1 7H?MR+01 
G 


3.0191 

0.0000 

3.H637 

0.0000 

2.3662 

0.0000 

1.7434 

0.0000 

1.4142 

0,0000 

1 .220R 

0.0000 

1.1034 

0.0000 

1.0353 

0.0000 

1.003R 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

01=  0.3H026PAR+00  0?=  0. 60F.3330E  +  00  03=  0. 1  ?5'J?0<)R  +  0? 

CnNVERGRNCt  (IR  LAGPiNGF  M(ILTPLIE« 

0.2ft7F>E+02  0.1A3HE  +  02  0.7190E+01  0.3S‘»5F  +  01  0.17Q7E+01  0.fl‘JR7E  +  00  0.4494E+00 
0.22A7F  +  00  0.n23E+00  O.HH^‘>E-Ol  O.^lllF-Ol 
FUNCTION 

0.‘>‘»S1F  +  01  n.SR3F.E+01  n.66fllE+01  0.S4fiSF  +  01  0.5133E  +  01  0.4f.F.3E4.01  0.3RR5E+01 
0.2SS4F+01  0.F>3M8F  +  00-O.HS0nE-01-O.U4  3i:-O2 

01=  n.3H02h24E  +  n0  02=  0.F)2F>4157E  +  00  03=  0. 1 6 1  25R 3E-t-n2 
CONVERGENCE  OF  LAGRANGE  MOLTPLiFR 
O.RlllE-01  0.1246E  +  00  O.mvF  +  OO  0.1322E+00 
FUNCTION 

-0,13F.0F+01-0.21  l‘»E  +  00-0. 72I2F-02-0.8722E-05 

01=  0.3R02624F+00  02=  0.62323900  +  00  03=  0. 1 3936030  +  02 
CONVERGENCE  OF  LAGRANGE  MULTPLIER 
0.13220+00  0.1433F+00  0.14410+00 
FUNCTION 

-0. 3209E +00-0. 1 6240-0 1-0.4292E-04 

01=  0.3H02624F+00  02=  0.61429600+00  03=  0.15360360+02 
CONVERGENCE  HE  LA&RANGF  MOLTPLIFR 
0.14410+00  O.13HO0+OO 
FUNCT  KIN 

0.1308E+00-0.42ROE-02 

01=  0.3H02624F+00  02=  0. 39H4214E  +  00  03=  0. 1473093F  +  02 
CONVERGENCE  UF  LAGRANG+  MOLTPlIFR 
0.13HOE+00  0.1263F+00  0.12690+00 
FUNCTION 

0.2R64F+00-0. 1629F-ni-0.44H2F-04 

01=  0.3H02624E+00  02=  0.5  7902370  +  00  03=  0. 1 525220E  +  02 
CONVERGENCE  UF  LAGRANGE  MOLTPLIFR 
0.12690+00  0.1330F+0n  0. 13330+00 
FUNCTION 

-0.2290E +00-0. « 3 400-02-0. 1  3  33F-04 

01=  0.3H02624E+00  02=  0.56296610+00  03=  0 . 1 6R o 37 R 0 +0 2 
CONVERGENCE  OF  LAGRANGE  miilTPLIER 
0.13330+00  0.13330+00  0.13760+00 
FUNCTION 

-O.3969F+OO-O.3123F-OI-O.43H70-O3 

01=  0.3H02624E+00  02=  0 . 554 16H9E  +  00  03=  0. 1  7  744440E +02 
CONVERGENCE  UF  LAGRANGE  MOLTPLIER 
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0.lb7f.F+00  0.1678fi  +  00  0.lftB3E  +  00 
HJMCl ION 

-0. 2  IS  14H  400-0. 10406-0 1-0.?  193e-0« 

01=  0.3H026?4fc-f00  02*  0.55004086+00  03=  0.18010446+0? 
C()NV68r.6NCE  06  LAGRANfip  MllLTPLlER 
0.16836+00  0.17116+00 
FDNCT ION 

-0.68266-01 -0.80116-03 

01=  0.38026746+00  02=  0.54810826+00  03=  0. 18070')lE+0? 
CONVERGENCE  1)6  LAGRANGE  MOLTRLIEH 
0.17116+00  0.17146+00 
FUNCTION 

-0.78456-02-0.17176-04 

01=  0.3802624E+00  02=  0.54716816+00  03=  0.18070446+02 
CONVERGENCE  06  LAGRANGE  MULTPLIER 
0.17146+00  0.1711E+00 
FUNCTION 

0.82306-02-0.16216-04 

01=  0.38026?4E+00  02=0.54669176+00  03=  0.18058246+02 
CONVERGENCE  06  LAGRANGE  MULTPLIER 
0.171 16+00  0.17076+00 
EUNCT ION 

0.9 362 6-0? -0.1 6? 16-04 


01=  0.38026?4E+00  02=  0.54644116+00  03=  0.18046256+0? 
CONVERGENCE  U6  LAGRANGE  MULTPLIER 
0.1707E+0C  0. 17046+00 
function 

0.7 164 6-02-0.76796-05 

01=  0.38026746  +  00  02=  0. 54630506+00  03=  0.18037486+02 
CONVERGENCE  06  LAGRANGE  MULTPLIER 
0. 1704E+0C 
FUNCTION 
0.48286-02 


NORM  SQUARE))  Ul 
O.OOOOOOOE+00 
0. 74442166+01 
;).H516630E+01 

RELATIVE  ERROR 
O.IOOOOOOE  +01 
0.232 166HE  +00 
0.19037721 +00 


0.0000000) +00 
0.68178121 +01 
0.681 76161 +01 


6 

0.71458726+01 

0.81453856+01 

0.85172266+01 


0.29244856+00 

0.20967556+00 

0.19035076+00 


0.68187076+01 

0.68175936+01 

0.6H1760HE+01 


0. 76003816+01 
0.84136126+01 
0.85173636+01 


0.26786146+00 
0. 19499066+00 
0.19034406+00 


0.68175636+01 

0.68180306+01 

0.68176136+01 


0.77801386+01 

0.84463516+01 

0.85173896+01 


0.25484826+00 

0.14120156+00 

0.14034196+00 


0.6817606E+01 
0. 68176116+01 
0.68175496+01 


GG=  0.5554)096+0?  66=  0.85114676+01  6=  0.19035956+00  0=  0 


6 

-0. 4309-, 98|  +00-0.45400566-01  0.87210796  +  00  0.44642066+00- 
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0.78924406+01 

0.85144756+01 

0.85115516+01 


0.24566516+00 

0.19050226+00 

0.14035456+00 


0.68218486+01 

0.68183886+01 

0.6812764E+01 

.68126986+01 


0.46832216+00 


-0.1137715E*01-0.1040341E+00  0.  144 14  1 '5F  +  01  0. 54Hfl ?4 1  E+00-0.  544344(16  +  00 
-0,45613556  +  00  0 .438S3 1  ?F  +  00  0.3?l‘l«?66+00-0.44?l  35«F+00-0.«>3fl556?6-01 
0.84307076+00-0. 33977576 +00-0. 75 104536+00  0.47939906+00  0.13424006  +  00 

FM 


0.290 

0.454  1 

.000  0.981  0.527  0.424  0.377 

0.548 

0.552  0.333 

FP 

-43.4 

157.5 

0.0  -135.5  84 

.7  -93. 

5  75.3 

-133.2 

14.1  144.0 

TF 

2.9092 

0.2108 

3.3419  -0.3124 

2.4455 

-1.1181 

1.0324 

-1.5195  -0.4792 

-0.9359 

-1.3893 

0.3440 

-0.8144  1.3489 

0.2045 

1.2814 

0. 4488 

0,5232  0.4525 

-0.0425 

0.1049 

-0.1550 

-0.0400  -0.0413 

-0.0735 

0.0098 

-0.1298 

-0.0157  -0.2001 

-0.0454 

-0.2308 

-0.1195 

-0.1982  -0.1489 

-0.1172 

-0.1224 

-0.0419 

-0.0587  -0.0272 

-0.0250 

-0.0701 

-0.0459 

-0.1221  -0.1474 

-0.1404 

-0.2071 

-0.1780 

-0.2410  -0.1303 

-0.2570 

-0.0185 

-0.1927 

0.0276  -0.0393 

-0.0583 

0.0494 

-0.1318 

-0.0044  -0.1094 

-0.0273 

-0.1478 

0.1114 

-0.4242  0.2525 

-0.4247 

0.2133 

-0.3445 

0.0923  0.5111 

0.1132 

1.7708 

0.2597 

TFM 

2.9148 

3.3544 

2.8905  1.8371 

1.1544 

1.4347 

1.5757 

1.2980  0.8491 

0.4547 

0.1872 

0.0575 

0.0742  0.1307 

0.2105 

0.2599 

0.2479 

0.1495  0.0721 

0.0370 

0.0942 

0.1914 

0.2419  0.2994 

0.2882 

0.1934 

0.0480 

0.0745  0.1319 

0.1129 

0.2014 

0.4937 

0.4420  0.3780 

0.5235 

1.7897 
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IX.  COMPUTATION  OF  SYNTHESIZED  PATTERNS  FROM  THE  SOURCE  VECTOR 


The  program  of  this  section  stores  the  synthesized  patterns  on  a 
direct  access  data  set  so  that  they  may  be  plotted  by  the  program  of 
section  X  of  this  report. 

The  activity  on  data  sets  1  (punched  card  input)  and  6  (direct  access 
input  and  output)  is  as  follows. 

RFADd.U)  N6,  N,  NX,  NP 

14  FORM.\T(20n) 

READ(1,16)(BK(I).  I-l,  NK) 

16  FORMAT (5E14. 7) 

READ(1,14)(NG(I),  1-1,  NK) 

READ(1,18)(X(I),  I-l,  N) 

18  FOR>UT(10F8.4) 

RFJ\D(1,18)(Y(1),  I-l,  N) 

REWIND  6 

IF(N6)  11,  11,  12 

12  DO  13  J-1,  N6 
READ  (6) 

13  CONTINUE 

11  FN  =  NP-1 

DO  25  JK-1,  NK 
KG  -  NG(,IK) 

DO  29  JG  =  1,  KG 
READ(1,16)(F(I),  I-l,  N) 

29  CONTINUE 

25  CONTINUE 

WRITE(6)(TFM(I),  1=1,  J9) 

The  X  and  y  coordinates  of  N  point  sources  are  read  in  as  X  and  Y.  A 

source  vector  f  at  propagation  constant  BK(JK)  is  read  in  through  F. 

Inside  DO  loop  29  the  synthesized  pattern  corresponding  to  the  source 

vector  F  is  calculated  as  the  magnitude  of  expression  (2-17)  evaluated  at 

2^ 

t  -  0,  All,  2A(Ji...  In  radians  where  Aii  =  — .  The  synthesized  patterns 
are  stored  one  after  the  other  in  TFM. 
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Mlnlnmm  allocations  are  given  by 


COMPLEX  T(N*NP),  F(N) 

DIMENSION  BK(NK).  NG(NK),  X(N) ,  Y(N} 

XCS(N*NP),  TFM(J9) 

where 

J9  -  NP*  [  NG(I) 

I-l 

Nested  DO  loops  23  and  24  put  cos  ^ -f  y^  sin  ((of  (2-17)  In  XCS. 

The  Index  JK  of  DO  loop  25  denotes  propagation  constant  BK(JK).  Nested 
DO  loops  26  and  27  store  the  matrix  [T]  of  (2-16)  by  rows  In  T  .  The 
Index  JG  of  DO  loop  29  denotes  the  JG^**  synthesized  pattern  to  be  computed 
at  propagation  constant  BK(JK).  Nested  DO  loops  31  and  32  store  the  syn¬ 
thesized  pattern  ||[T]^||  of  (2-17)  In  TFM. 

The  sample  punched  card  Input  data  Is  such  that  the  patterns  syn¬ 
thesized  by  the  programs  of  sections  II,  III,  V,  VI,  VII,  and  VIII,  of 
this  report  are  evaluated  at  every  2-1/2  degrees  and  stored  on  Che  first 
record  of  data  set  6. 


63 


LISMNG  OF  PROGRAM  TO  rOMPOTF  SVNTHESI^Ff)  PATTERNS  FROM  SOURCE  VECTOR  E 


//  (OOTA,EEtlSS,n,*MAI(TZtJnE'  ,  R  EG  I  nN«?OOK 

//  FXFC  WATFIV 

//GO.FTOhFnni  1)0  r)SNAMF  =  FE0034.REVl  fniSPBnLf)fllNIT*33T0,  X 

//  vnLUMF=SER=Stl0009,nCR*(RECFM*VS,RLKSI7E«?5‘)A,LRECL  =  ?5P?tX 

//  HIJFNO*!  ) 

//GO. SYS  IN  00  * 

JJflh  MAIJTZ  ,TlME  =  ltEAr,ES*40 

COMPLEX  0,01 , T (  14S0I  tFl 10( 

DIMENSION  HK ( 1 01  ,NG(  10> .XIIO) .V( 10) *XCS( 1450 )  ,TFMI 2P00) 

0= ( O. . I  .  I 
PI=  1.141593 
RFAI)(1,14)  N6,N,NK,NP 

14  F0R^AT(?0I3) 

WR  I  fEl 3, 15 )  N6,N,NK,NP 

15  FOR^All'O  N6  N  NK  NP‘/1X,4I3) 

REA  )(  1  .1  />)  (  BK(  I  )  .  I  =  1.NK) 

15  FnR«1AT(5E14.7  I 

WR I  IE  (  3,17 )  (BKI I  )  .I«1«NK) 

17  FORMAT! 'OBK'/( 1X,5E14.7)) 

REA  )( I • 14) (NG( I  )  ,  I<1,NK) 

WR I rE(3,2fl )  (NG(I  )  ,I  =  1,NK) 

2H  FORMAT! 'ONG '/( 1 X, 2013) ) 

READ! 1 , IR ) IX !  I  )  ,  I  :)  ,N) 

M  FORMAT!  lOf  fl.4) 

WR1TF(3,19) IX!  I  )  ,I=1,N) 

19  FORMAT!  '0X'/!lX,l0Ffl.4)) 

RFA0(1,18)IY!I  ).I=1,N) 

WR ITfc (3,20) !Y !  I) .  I  =  1,N) 

20  FORMAT! 'OY '/ I  IX, inFfi.4) ) 

REWIND  5 

IFIN6)  11,11,12 

12  00  13  J*1,N6 
RFAD(6I 

13  COATINOE 
11  FN=NR-1 

0FL=2.oPI/FN 
J1  =  0 

00  23  ,)  =  1,NP 
Sl= !J-1 IfOEL 
CS=COS! SI ) 

SN=SIN!S1 I 
no  24  1=1, N 
Jl: Jl+1 

xc;  !J1  l  =  X!  I  )4CS*Y|I  )«SN 

24  COOTINUfc 
23  COfiTINUF 

,|9;rO 

DO  25  JK=1,NK 
KG  iNGIJK ) 

)3H  :RK  (  JK  ) 

Jl  =0 

on  25  J=1,NP 

DO  27  I  =  1  ,N 

Jl  :J1  +  1 

S  1  M3HYXCS ! Jl  I 

I  !  )1 )=COS! SI )+0*SINISl  ) 

2  /  (.(HU  NOE 

25  CO  UINOE 


nil  2»  jr.=i,KO 

READ! 1 |16| (F(  I  l,l  =  l,N) 

MRITE(3,30IIF(l)«l<ltN) 

30  F0RMAT('0F'/(lX,iHlA.7)  1 
J1»0 

nn  31  .i»iiNF 

1111=0. 

on  32  Iti.n 

Ji=Ji+i 

lll*Ul+t(Jl  )*F(  I  ) 

32  CDNTINIIF 
JR'JR+l 

TFM(  JR  )=CAHS  llJl  ) 

31  CONTINUE 
29  CONTINUE 
25  CONTINUE 

WR1TE(6)ITFM( I ),I=1,J91 

N9=J9/NP 

Jl*l 

DO  33  J=1,N9 
J2*Jl+9 
WRITF(3,35)  J 

35  FORMAT! 'OFIRST  10  POINTS  ON  THE*. 13.'  TH  PATTERN') 
WRITE(3,34) (TFM(  I  I, r=Jl  , J2) 

34  F0RMAT(1X,10FH.4) 

Jl-Jl+NP 

33  CONTINUE 
STOP 
END 


SOATA 

0  10  1145 

0.290flflfl2E+00 
4 

0.0000  5.340H  10.4048  15.5077  19.2448  19.7448  15.5077  10.4048  5.3408  0.0000 
-lO.OOOn  -9.6320  -8.4706  -6.2894  -7.5559  7.5559  4.2894  8.4704  9.4320  10.0000 
-0.5421209E+00  0.5724890E-01  0.  1 187R10E  +  01  0. 5949R31E+00-0. 1447048E  +  01 
-0.944R900E  +  00  0 .6557452F  +  00  0 . 14  37544e+01-0. 3453 1 758 *00-0. 7954181E+00 
-0.9229400E-02  0. 1055890F*01-0. 3044195E  +  00-().  R304R94E  +  00  0. R837318E-01 
0. 12472 18E+01-0. 738 1020F  +  nO-0.8539429F+00  0.5231934F+00  0.7333544E-01 
-0.4131346E  +  00-0.2757374E+00  0.  1744744E+ni  0. 9942800E  +  00-0. 1 532R91E  +  01 
-0.1988470E  +  01  0 .47873 1 8E  +  nO  0 . 2507443E+01  0.?8 1 1 250E +00-0 . 1595 194E+01 
-0.3373244E+00  0. 1369034E  +  01  0. 223438 lE+00-0.  14355R2E  +  01  0. 1 407297E-01 
0.1491R13E+01-0.4139780E+00-0.1279405E+01  0.4B12400E+00  0, 3544 1 19E+00 
-0.3427129E+00  0 . 2398499E +00  0. 7955879E+00  0.  1R84414E  +  00-0.9377459E+00 
-0.3302 187E+00  0 . 1471 369F  +  00  0 .7894499E+00-0.4094427E-03-0, 7  1  24  lRRE+00 
-0.3034758E+00  0. 449  7439F +00-0. 1 449025E-01-0.  1 821 197E  +  00-0. 1 9431 34F+00 
0.41 73351 E+00-0. 5551 674E  +  00-0.4373773E+00  0.428 1  34BE +00-0. 1 71 7409E+00 
-0.4177318E+00-0.4299459E-01  0. 844 1 843E+00  0. 4310744E  +  n0-0. 94 1  2030F+00 
-0.1100197E  +  01-0.9487757E-01  0 . 14  14402E+01  0. 5494794E +00-0. 52932R4F+00 
-0.440R940E+00  0 . 474740PF +00  0. 3 1 OH431E+00-0.  448R897E  +  00-0. 8R91094E-01 
0.81 72523E  +  00-0. 33099  10F  +  00-0.7249035F+00  0. 444701 8F +00  0. 1289077E +00 
-0.3523402E+00  0. 246378 2F +00  0. 8 1  79 1 lOF+00  0. 1 94 1 044E  +  00-0. 9440R98E  +  00 
-0.3397281E+00  0 . 1 669434F+00  0.81 1 7170F+00-0.4355443E-03-0.7188425E+00 
-0.31 1942 lE+00  0.4825977F+00-0.1777355F-01-0. 1 R75400F+n0-0. 701 4137F+00 
0. 4348063E+00-0. 5704 170E +  00-0. 444  594 lE+00  0.44n0374E +00-0. 1743947E+00 
-0.4309598E  00-0 . 4540054E-01  0.8721079E  00  0.4449204F  00-0.9483771E  00 
-0.1137715E  01-0.1040341F  00  0.1441415F  01  0.548H741E  00-0 . 5443448E  00 
-0.4541355E  00  0.43893178  00  0.3219824E  00-0. 4621 35HE  00-0.93R5547F-01 
O.H430707E  00-0.33977578  00-0 .75 10453F  00  0.4793990F  00  0.1342400E  00 


SSTOP 
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/« 

// 


MH  FNTI-D  UU1PUI 

Nh  N  NK  NP 
(1  10  1146 


HK 

0.2<J0HHH2P  +  0f' 


X 

0.0000  6.^(.0H  10.604H  15.50P7  l‘J.?A4ft  14.244R  15.50?7  10.7.04fl  5.340R 

Y 

-10.0000  -9.6)20  -fl.4706  -6.2H94  -2.5569  2.5559  6.2R94  R.4706  9.6320  1 

P 

-0.5421209t+0(i  0.5224fl90P-01  0. 1 1 H7fl lOF+01  0. 5969R3 1 E  +  00-0. 1 46704RP  +  01 
-0.944H9n0E  +  OO  0.65524  52P  +  00  0. 14  32544F  +  01-0. 36531  75P+00-0. 7956 IR lE+00 
-0.9229600P-0:'  0.1056fl90F+01-0.3044195E+00-O.R306H94F  +  00  0. flH373 1  RE-01 
0.1247218F+0  -0.7381020E+00-0.8539429E+00  0.5231936E+00  0.2333546E-01 


-0.6131346L-40(I-0.2752376F+00  0. 1244744E+01  0.9962R00E+OO-0.  1  532R91E+01 
-0.19H8470F  +  0  0.47873)8F  +  O0  0 .25O2643F4'01  0 . 2R 1 1 250F +00-0, 1 595 194E +0 1 

-0.3373244F  +  00  0 . 1 369034F +0 1  0.22363R1E+00-0. 1435582E  +  01  0. 1 607297F-0 1 
0.1691H13E+0  -0.4139780F+00-0.1279405E+01  0.4812400E+00  0. 3546 1 19E +00 


-0.3427129P+0()  0 . 239H499E +00  0. 7955879E+00  0. 1 RR641 6F+00-0. 9377459E+00 
-0.33021 8 7e  +  00  0 . 1  62 1 369F  +  00  0.7R94699E+00-0.4096627E-03-0.21 261RRE+00 
-0.303675HE+00  0 . 4692439E +00-0. 1669025E-01 -0. 1 R21  197E  +  00-0. 1 963 1 36E+00 
().6l7  3351E  +  Oi)-0.5551624E  +  00-0.4323723F+00  0 .42R 1 34RE +00-0 . 1  7 17609E+00 


F 

-0.4177318F+0()-0.4299659E-01  O.H461R63E+00  0. 431 0746E +00-0. 941 2030F+00 
-0.1100197E  +  0 ,-0.9  6H77  62E-01  0 . 14 16602F+01  0 . 5496 796F +00-0. 52932R6E +00 
-0.440H940b  +  00  0.4267  602E+00  0. 3 10863 lE  +  00-0. 44fl8R97F +  00-0. 8  891 094E-01 
O.H172523E+n )-0.3309910E+00-0.7269035E+00  0.464701RF+00  0. 12R9072E+00 


F 

-0.  )523602E  +  0i)  0.2463782E+00  0.81791 lOF+00  0. 1941064F  +  00-0,9640R9RE+00 
-0. 3  39728 lE  +  O)  0 . 1 669634E  +  00  0 .8 1 1  71 20F+00-0.6355643E -03-0 . 2 1 RR625E+00 
-O.3119421E  +  00  0.4825922F+00-0.1727355F-01-0.  1R75600E  +  00-0.2016132F+00 
0. 634806 3E  +  0 1-0.57061  20F  +  00-0, 444596 lF+00  0. 44003 74E +00-0. 1 763967E +00 


F 

-().430959RE  00-0.4540056E-01  0.8721079F  00  0.4469206F  00-0. 968322  IE  00 
-0.1137715E  01-0.1040341E  00  0.1461415E  01  0.56R8261E  00-0.544344RE  00 
-().4561355E  0)  0.4389312F  00  0.3219R26E  00-0. 4621 35RF  00-0. 93R5562F-01 
0.8430707E  03-0.3397757F  00-0 .75 10453E  00  0.4793990E  00  0.1342400E  00 

FIRST  10  POINI.S  ON  THE  1  TH  PATTERN 

1.9967  2.2T96  2.5691  2.7961  2.9730  3.0940  3.1553  3.1555  3.0953 

FIRST  10  PfllNfS  ON  THE  2  TH  PATTERN 
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i 


O.OOOf 

0.000( 


2.977 


2.7330  3.0170 

3.2373  3.4001 

3.4942 

3.5161 

3.4653 

3.3446 

3.1597 

FIKST  10 

POINTS  ON  THE  3 

TH  PATTERN 

1  .«»3«J7 

2.190H  2.4135 

2.6013  2.7488 

2.8520 

2.9081 

2.9160 

2.8760 

2.7901 

FIRST  10 

POINTS  ON  THE  4 

TH  PATTERN 

2.3337 

2.5976  2.8257 

3.0106  3.1459 

3.2272 

3.2516 

3.2186 

3.1295 

2.9881 

FIRST  10 

POINTS  ON  THE  5 

TH  PATTERN 

1 .9R35 

2.2516  2.4804 

2.6734  2.8250 

2.9311 

2.9888 

2.9969 

2.9558 

2.8674 

FIRST  10 

POINTS  ON  THE  6 

TH  PATTERN 

2.A0P9 

2.6813  2.9168 

3.1076  3.2473 

3.3312 

3.3564 

3.3223 

3.2304 

3.0844 
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X.  PLOTS  OF  SPECIFIED  AND  SYNTHESIZED  PATTERNS 


The  program  of  this  section  plots  the  synthesized  patterns  stored 
on  direct  access  data  set  6  by  the  program  cf  section  IX  of  this  report. 

Tie  activity  on  data  sets  1  (punched  card  Input)  and  6  (direct 
access  Input  and  output)  Is  as  follows. 


READ(1,10)  N6,  NG,  NP,  M, 

SCL 

10 

FORMAT(4I3,  E14.7) 

READ(1,12),  (PHI(l),  1-1, 

M) 

12 

FORMAT (10F8, 4) 

READ(1,12)(H(I),  1=1,  M) 

J9  =  NP*NG 

REWIND  6 

IF(N6)  16,  16,  17 

17  no  18  J-1,  N6 
READ (6) 

18  CONTINUE 

16  READ(6) (TFM(I),  I«l,  J9) 

The  specified  pattern  evaluated  at  M  angles  PHI  In  degrees  Is  read  into  H. 
NG  synthesized  patterns  evaluated  at  NP  equally  spaced  angles  ip  are  read 
into  TIM.  All  patterns  are  multiplied  by  the  scale  factor  SCL  and  then 
plottec  In  Inches. 

Minimum  allocations  are  given  by 

DIMENSION  PHI (M)  ,H(M)  ,XM(M) , YM(M) , 

TFM(NP*NG),  CS(NP),  SN(NP),  X(NP),  Y(NP) 

DO  loop  15  stores  the  x  and  y  coordinates  of  the  specified  pattern 
in  XM  and  YM  respectively.  DO  loop  21  stores  in  CS  and  SN  the  cosine  and 
sine  of  the  angles  at  which  the  synthesized  patterns  are  evaluated. 

"he  Index  JG  of  DO  loop  20  denotes  the  JG^^  synthesized  pattern 
to  be  >lotted.  DO  loop  22  puts  tick  marks  on  the  x  axis  drawn  by  statement 
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26.  DO  loop  23  puts  tick  marks  on  the  y  axis  drawn  by  statement  27.  DO 
loop  24  plots  the  first  10  points  on  the  specified  pattern.  The  first  10 
points  are  sufficient  because  the  particular  specified  pattern  used  Is 
xero  at  all  but  the  first  9  points.  For  a  general  specified  pattern,  DO 
loop  24  should  be  changed  to  run  from  J  •  1  to  J  ■  M.  Statement  28  draws 
straight  lines  between  the  x  and  y  coordinates  of  the  synthesized  pattern 
prepared  by  DO  loop  25. 
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LISTING  OF  PROGRAM  TO  PLOT  SPECIFIED  AND  SYNTHESIZED  PATTERNS 


//  (0034,EE.30S»l. t 12l,*MA(ITZtjnE •,ReGI0N=140K 

//  MSG  T,  70  INCHES  OF  PLOT  PAPER  IS  REOOIRED 
n  EXEC  FOKTGCLG 
//FORT.  SYS  IN  DO  * 

DIMENSION  PHI  (  36  I  ,  H  (  36  )  ,  XM  (  3f>  >  ,  YM I  3N )  »  TFM I  ?900  I .  XX  (  4  )  ,  YY  ( 4  ) 
DIMENSION  CS  (  14S  )  ,  SN  (  14S  )  .  ARF  A  (  400  I  ,xn4S),Y(14St 
PI=3.1415<J3 
C=P| /IRO. 

READ (1,101  N6,NG,NP,M,SCL 

10  FrRMAT(4I3,ei4.7) 

WBITE(3,11I  N6,NG,NP,M,SCL 

11  FfRMATCO  N6  NG  NP  M  •  ,5X  ,  •  SCL  •  /  1 X  ,41  3  ,E  14 . 7  ! 

READ! 1 , 12) (PHI (I ),I=1,M) 

12  FORMA! ( 10FH.4) 

W6  I  TK  3, 13)  (PHI  (  I  I  ,  I=1,M> 

13  FORMAT (  'OPHI  •/( 1 X, 10F8.4  )  ) 

Rf A0( 1 ,12) (H( I ),I=1,M) 

WMTE(3,14)(H(1).I=1,M) 

14  FT RMAT (  'OH'/ ( 1X,10FH.4  )  ) 

DO  IS  J=1 ,M 

S1=C*PHI (J) 

SJ=SCL*H(.)) 

XR ( J)=6.+S2»CnS(Sl ) 

Yf)l  J  )=S.  +  S?*SIN(  SI  ) 

15  Cl  NTINIIF 
J‘'  =  NP>i‘NG 
REWIND  6 

!) (N6)  16,16,17 
17  Of  la  .1=1, N6 
Rl A0(6  ) 
la  cr  NTINOE 

16  K)  AD(6)  (TFM(  I  ),I=1  ,.)9) 

Wt  ITF(3.1R)(TFM(I  ),I=l,in) 

1<J  FORMAT  (  'OTFM '/(  1  X,  lOFa.4  )  ) 

X)  (  1  )  =  1  . 

X  •;  ( 2  )  =  1  1  . 

Y' (11=5. 

Y  '(2 )=5. 
x;  (3)  =  6. 

X  ((4 )=6. 

Y''(3)=0. 

YM4)  =  10. 

Fri  =  NP-l 
D ■L=2 .*PI /FN 
00  21  J--1,NP 
S  l  =  (.)-l  I'i'DEL 
C;,(  J)=C0S(S1  ) 

S0( J I=SIN(S1 ) 

21  CONTINUE 
OLL  PLOTID 
C.tLL  PLOTS!  AREA, 400  I 
J  (  =  0 

DO  20  JG=1,NG 

26C\LL  LINE(XX(1 ),YY(1 ),?,!, 0,0) 

S ,  =  11 . 

m  22  .)  =  !  ,11 

f.  \LL  SYMB()L(S1 ,5.  ,.14,13,0.,-ll 
S1  =  S1-1  . 
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22  CDNTINllfc 

27  CALL  LINe(XX(3)tYy(3(,?,liO.OI 
SI  *10. 

nn  23  .1=1,11 

CALL  SYMHnL(h.,Sl,.14,13,'»0.,-l) 

S1=S1-1  . 

23  CONTINOt 

DO  24  .1  =  1 ,10 

CALL  SYMROL(XM(.)(,YM(J),.  14,4,0.,-!  » 

24  CONTINUE 

on  25  .1=1, NP 
.)1=J1  +  1 

S!=SCL*TFM( J1  I 
X(J)=6.  +  S1*CSI.I) 

Y(  J)  =  5.  +  S1*SN(.I) 

25  CONTINUE 

28  CALL  LINe(X,Y,NE, 1,0,0) 

CALL  PL0T(ll.,0.,-3) 

20  CONTINUE 

CALL  PLnT(3.,0.,-8) 

STOP 

ENO 

/♦ 

//r,0.ETn6E001  on  nSNAME=EE0034.REVl  ,nTSP=nLD,UNlT  =  3330,  X 

//  VOLUME =SEK =500009, OCR=(RECFM=VS,RLKS I 7E=25Q4,L«ECL=?54?,Y 

//  RUFNU=1 ) 

//GO. SYS  IN  on  = 

0  6145  36  0.5000000F+00 

5,0000  15.0000  25.0000  35.0000  45.0000  55.0000  65.0000  75,0000  85.0000  95.0000 
105. 0000115.000012  5.0000135.0000145.0000155.0000165.0000175.0000185. 0000 195. 0000 
205.0000215.0000225.0000235.0000245.0000255.0000265.0000275.0000285.0000295,0000 
305.0000315,0000325.0000335.0000345.0000355.0000 


5,0191 

3.8637 

2.3662 

1.7634 

1.4142 

1.2208 

1.1034 

1.0353 

1.0038 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

/» 

// 

PRINTED  OUTPUT 

N6  N(;  NP  M  SCL 

0  6145  36  0.5000000E'»00 

PHI 

5.0000  15.0000  25.0000  35.0000  45.0000  55.0000  65.0000  75.0000  85.0000  95.0000 
105.0000115.0000125.0000135.0000145.0000155.0000165.0000175.0000185.0000195.0000 
205.0000215.0000225.0000235.0000245.0000255.0000265,0000275.0000285.0000295.0000 
305.0000315.0000325.0000335.0000345.0000355.0000 


H 


5.0191 

0.0000 

0.0000 

0.0000 

3.H637 

0.0000 

0.0000 

0.0000 

2.3662 

0.0000 

0.0000 

0.0000 

1.7434 

0.0000 

0.0000 

0.0000 

1 .4142 
0.0000 
0.0000 
0.0000 

1.2208 

0.0000 

0.0000 

0.0000 

1.1036 

0.0000 

0.0000 

1 .0353 
0.0000 
0.0000 

1.0038 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

TFM 

1 .9967 

2.2996 

2.5691 

2.7961 

2.9730 

3.0960 

3.1553 

3.1555 

3.0953 

2.9777 
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